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The small- a; deep inelastic scattering in the saturation region is governed by the non-linear evo- 
lution of Wilson-line operators. In the leading logarithmic approximation it is given by the BK 
equation for the evolution of color dipoles. In the next-to-leading order the BK equation gets con- 
tributions from quark and gluon loops as well as from the tree giuon diagrams with quadratic and 
cubic nonlinearities. We calculate the gluon contribution to small-x evolution of Wilson lines (the 
quark part was obtained earlier). 
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A general feature of high-energy scattering is that a fast particle moves along its straight-line classical trajectory 
and the only quantum effect is the eikonal phase factor acquired along this propagation path. In QCD, for the fast 
^ quark or gluon scattering off some target, this eikonal phase factor is a Wilson line - the infinite gauge link ordered 

(~| along the straight line coUinear to particle's velocity n^: 
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I. INTRODUCTION 



U'^{x±) ^PexpUg f du Tif, A^'{un + xj_)\, (1) 

J — oo 



Here is the gluon field of the target, x± is the transverse position of the particle which remains unchanged 
throughout the collision, and the index i] labels the rapidity of the particle. Repeating the above argument for the 
target (moving fast in the spectator's frame) we see that particles with very different rapidities perceive each other as 
Wilson lines and therefore these Wilson-line operators form the convenient effective degrees of freedom in high-energy 
^y-^ QCD (for a review, see ref. fV). 

Let us consider the deep inelastic scattering from a hadron at small xb = / {2p-q). The virtual photon decomposes 
into a pair of fast quarks moving along straight lines separated by some transverse distance. The propagation of this 
quark-antiquark pair reduces to the "propagator of the color dipole" U{x±)W {y±) - two Wilson lines ordered along 
the direction coUinear to quarks' velocity. The structure function of a hadron is proportional to a matrix element of 
f — . this color dipole operator 

9 . 1 . 

> U^x^,y^)^l-—Tr{U^^{x^)U^^^{yj^)} (2) 

switched between the target's states {Nc = 3 for QCD). The gluon parton density is approximately 

xbG(xb,m' = Q') ^ {p\U'^{x^,0)\p) (3) 

where rj = In . (As usual, we denote operators by "hat" ) . The energy dependence of the structure function is 
translated then into the dependence of the color dipole on the slope of the Wilson lines determined by the rapidity 77. 

Thus, the small-x behavior of the structure functions is governed by the rapidity evolution of color dipoles [21 [3]. 
At relatively high energies and for sufficiently small dipoles we can use the leading logarithmic approximation (LLA) 
where Ofg ^ 1, as Inxs ^ 1 and get the non- linear BK evolution equation for the color dipoles [HIS]: 



^^U{x,y) = J d z _ _ ^^^ [U{x, z) +U{y, z) - U{x,y) - U{x, z)U{z,y)] (4) 
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The first three terms correspond to the linear BFKL evolution [5] and describe the parton emission while the last 
term is responsible for the parton annihilation. For sufficiently high xb the parton emission balances the parton 
annihilation so the partons reach the state of saturation [7] with the characteristic transverse momentum Qg growing 
with energy 1/xb (for a review, see f5^) 

As usual, to get the region of application of the leading-order evolution equation one needs to find the next-to- 
leading order (NLO) corrections. In the case of the small-x evolution equation Q there is another reason why NLO 
corrections are important. Unlike the DGLAP evolution, the argument of the coupling constant in Eq. Q is left 
undetermined in the LLA, and usually it is set by hand to be Qs- Careful analysis of this argument is very important 
from both theoretical and experimental points of view. From the theoretical viewpoint, we need to know whether the 
coupling constant is determined by the size of the original dipole |a; — ?/| or of the size of the produced dipoles |a; — z| 
and/or |z — since we may get a very different behavior of the solutions of the equation On the experimental 
side, the cross section is proportional to some power of the coupling constant so the argument determines how big (or 
how small) is the cross section. The typical argument of as is the characteristic transverse momenta of the process. 
For high enough energies, they are of order of the saturation scale Qs which is ~ 2 3 GcV for the LHC collider, so 
even the difference between a{Qs) and a(2Qs) can make a substantial impact on the cross section. The precise form 
of the argument of as should come from the solution of the BK equation with the running coupling constant, and the 
starting point of the analysis of the argument of as in Eq. Q is the calculation of the NLO evolution. 

Let us present our result for the NLO evolution of the color dipole (hereafter we use notations X = x—z, X' = x — z', 
Y = y — z, and Y' = y — z') 
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Here (x is the normalization point in the MS scheme and h — — |n/ is the first coefficient of the /3-function. 

The result of this paper is the gluon part of the evolution, the quark part of Eq. (jsj) proportional to nj was found 
earlier [9j[l0]. Also, the terms with cubic nonlinearities were previously found in the large- A'c approximation in Ref. 
The NLO kernel is a sum of the running-coupling part (proportional to b), the non-conformal double-log term 

~ In ^^^Zs'zY [^-t)^ ^^'^ t^*^ three conformal terms which depend on the two four-point conformal ratios x''^y'^ ^^'^ 

'•^~^2y/lj^ "* ■ Note that the logarithm of the second conformal ratio In ^'^~^2y<I^ is absent. 

It should be emphasized that the NLO result itself does not lead automatically to the argument of coupling constant 
as in Eq. |4] In order to get this argument one can use the renormalon-based approach [12j: first get the quark part 
of the running coupling constant coming from the bubble chain of quark loops and then make a conjecture that the 
gluon part of the /3-function will follow that pattern. The Eq. ^ proves this conjecture in the first nontrivial order: 
the quark part of the /? - function |rij calculated earlier gets promoted to full b. The analysis of the argument of 
the coupling constant was performed in Refs. [HI [10] and we briefly review it in Sect. 7 for completeness. Roughly 
speaking, the argument of as is determined by the size of the smallest dipole min(\x — y\,\x — z\,\y — z\). 

The paper is organized as follows. In Sect. II we remind the derivation of the BK equation in the leading order in 
as- In Sect. Ill and IV, which are central to the paper, we calculate the gluon contribution to the NLO kernel of the 
small- a; evolution of color dipoles: in Sect. Ill we calculate the part of the NLO kernel corresponding to one-to-three 
dipoles transition and in Sect. IV the one-to-two dipoles part. In Sect. V we assemble the NLO BK kernel and in 
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FIG. 1: Leading-order diagrams for the small- a; evolution of color dipole. Gauge links are denoted by dotted lines. 



Sect. VI we compare the forward NLO BK kernel to the NLO BFKL results [T3]. The results of the analysis of the 
argument of coupling constant are briefly reviewed in Sect. VII. Appendix A is devoted to the calculation of the 
UV-divergent part of the one-to-three dipole kernel and in Appendix B we discuss the dependence of the NLO kernel 
on the cutoff in the longitudinal momenta. 



II. DERIVATION OF THE BK EQUATION 



Before discussing the small-x evolution of color dipole in the next-to-leading approximation it is instructive to recall 
the derivation of the leading-order (BK) evolution equation. As discussed in the Introduction, the dependence of the 
structure functions on xb comes from the dependence of Wilson- line operators 



U''''{xi_) = Pexp|i5 / du A^{u'n + a;_L)|, 



(6) 



on the slope of the supporting line. The momenta pi and p2 are the light- like vectors such that q = Pi — xbP2 and 

2 

P = P2 + where p is the momentum of the target and m is the mass. Throughout the paper, we use the Sudakov 
variables p = api + (3p2 + p_l and the notations = a;^^^ and = x^P2 related to the light-cone coordinates: 

3j i); — S /2 , X% — X S j 2 . 

To find the evolution of the color dipole (|2]) with respect to the slope of the Wilson lines in the leading log 
approximation we consider the matrix element of the color dipole between (arbitrary) target states and integrate 
over the gluons with rapidities rji > r] > r]2 = iji — A?7 leaving the gluons with r; < 772 as a background field (to be 
integrated over later). In the frame of gluons with 1] rji the fields with 77 < 772 shrink to a pancake and we obtain 
the four diagrams shown in Fig. [T| Technically, to find the kernel in the leading-ordrer approximation we write down 
the general form of the operator equation for the evolution of the color dipole 

^Tt{U,UI} = XLoTr{C/,i7t} + ... (7) 

(where dots stand for the higher orders of the expansion) and calculate the l.h.s. of Eq. ([t]) in the shock-wave 
background 
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In what follows we replace 



Shockwave 



by 



for brevity. With future NLO computation in view, we will perform 



the leading-order calculation in the lightcone gauge p-^Ap, = 0. The gluon propagator in a shock- wave external field 
has the formflU [H] 



{Al{x)Al{y)) = 9{x^y.)S- 
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FIG. 2: Leading-order diagrams proportional to the original dipole. 



where 



df,Ak) = gi, - —{k^p2„ + k^P2t,) ~ —P2t,P2u (10) 
^ sa sa 

Hereafter we use Schwinger's notations {x±\F{p±_)\y±^ = /(fp e'-'^P^'^^y'^^ F{pj^) (the scalar product of the four- 
dimensional vectors in our notations is x ■ y — -{x^,y, + x^,y,) — {x,y)±). Note that the interaction with the shock 
wave does not change the a-component of the gluon momentum. 
We obtain 

Formally, the integral over a diverges at the lower limit, but since we integrate over the rapidities ?7 > 772 we get (in 
the LLA) 

f° dv {A:{un + x^)Ai{vn + y^))F,^m. = - ia^Ar^ix^l^U'''' ^\y^) (12) 
Jo J-00 ^ PI PI 

and therefore 



The contribution of the diagram in Fig. [ijD is obtained from Eq. ( 13 1 by the replacement f^Ux ® t^U^ — > 

X y and the two remaining diagrams are obtained from Eq. I^Dy taking y — x (Fig. [ij;) and x — y (Fig. [l|i). 

Finally, one obtains 



" 7 (a;-z)^ TT^ "J (y-2)l 

so 

(Tr{^7.{7t})^l^[^ = ^/ ^'^^ _%2jy't [Tr{[/.^l}Tr{£/.[/t} - i-Tr{C/.t/t}]"^ (15) 

There are also contributions coming from the diagrams shown in Fig. [2] (plus graphs obtained by reflection with 
respect to the shock wave). These diagrams are proportional to the original dipole Tr {UxU^} and therefore the 
corresponding term can be derived from the contribution of Fig. 1 graphs using the requirement that the r.h.s. of 
the evolution equation should vanish for x — y since lim^^^y '§^"^^{^x^y} ~ easy to see that this requirement 

leads to 

{Tr{UxUl}r = ^ / d'z± -%i{y't z)l [Tr{t/.£/l}Tr{[/,t/t} - jV^Trlt/.t/t}]"- (16) 
which is equivalent to the BK equation for the evolution of the color dipole Q. 
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FIG. 3: "Cut self-energy" diagram. 



III. DIAGRAMS WITH TWO GLUON-SHOCKWAVE INTERSECTIONS 



A. "Cut self-energy" diagrams 



In the next-to- leading order there are three types of diagrams. Diagrams of the first type have two intersections 
of the emitted gluons with the shock wave, diagrams of the second type have one intersection, and finally there are 
diagrams of the third type without intersections. In principle, there could have been contributions coming from the 
gluon loop which lies entirely in the shock wave, but we will demonstrate below that such terms are absent (see the 
discussion at the end of Sect. VI I. 



For the NLO calculation we use the lightcone gauge P2^^ — 0. Also, we find it convenient to change the prescription 
for the cutoff in the longitudinal direction. We consider the light-like dipoles (in the pi direction) and impose the 

cutoff on the maximal a emitted by any gluon from the Wilson lines so 



= Pexp 



dup'^A^iupi + XA_) 



Al{x) = /cf4fc0(e''-|afe|)e-*-A^(fc) 



(17) 

As we will see below, the (almost) conformal result ([s]) comes from the regularization (17 1. In Appendix B we will 
present the NLO kernel for the cutoff with the slope ([6]). We start with the calculation of the Fig. [3^ diagram. 
Multiplying two propagators ([9]), two 3-gluon vertices and two bare propagators we obtain 



g2 y du \ dt>(A^(upi + Xi_)Al{vpi + y_L)) 
1 as' 



(18) 



4ai(a - ai);7r'f/15'e-'(«i-'=i'^)^-'(9^-'=^'^')^ d,x{api + (3p2 + qii_ + ki^) dy,{api + (3'p2 + q2± + fc2±) 



iP-pi-p2 + ie){f3' -l3[-p'^ + ie){(i - ie){P' - le) a(3s - (gi + 92)! + it 



a(3's - {ki + k2)\ + ie 



d^dc'iPi + P1P2 + qi±) d^AaiPi + P1P2 + ^I-l) d^^{{a - ai)pi + (32P2 + q2±) d'^^,{{a - ai)pi + /3^p2 + fc2±) 



ai/3is - qfj_ 



ai(3[s — fc^j_ + ie (a ~ ai)P2S — (7|j_ + it 



(a — ai)f3'2S — fc|_|_ + ie 



^''^{api + qi±, {a - ai)pi + q2±,-api - qij_ - q2±) T^"" ^ (api + ki±, {a - ai)pi + k2±,-api - kij_ - k2±) 
where 

r^i.A(p, k, ~p~k)^{p- k)xg^u + (2k + p)f_,g^\ + {-2p - k)i,gx^_, (19) 
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In this formula ^g^^ comes from the integration over u parameter in the l.h.s. and from the integration of 

the left three-gluon vertex over the half-space a;* > 0. Similarly, we get from the integration over v parameter 

and 757 — ai^ Ri j_- from the integration of the right three-gluon vertex over the half-space cc* < 0. The factor 1 in the 
r.h.s. is combinatorial. Note that in the light-cone gauge one can always neglect the (3p2^ components of the momenta 
in the three-gluon vertex since they are always multiplied by some d^^j. 
Taking residues at /? = /9' = and P2 = —Pi, P2 = ~Pi we obtain 



(20) 



^g^jT^'f"''' Jdadaid(3id(3[ J d^zd^z' J d^qid^q2d^hd^k2 e*(9i+92-)^-»(fci+fe2,y)x 



aija - ai) ^jn,i' rrll'^i(qi~ki)z^i(q2~k2)z' (glJ- + g2±)A (fci± + k2±)\' 

a2 ^-'^ (<7i + g2)i (fcl + fc2)i 

d^{aipi+ qij_) d^^,{aipi + ki±) d,,'' {{a - ai)pi + q2±) - ai)pi + fc2±) 



ai/3is — qfj_ + it ail3'iS — k\_^ + ie —(a — ai)(3is — q^j^ + it —{a — ai)(3iS — -|- it 

V^''^{aipi + (a - ai)pi + 172^, -api - qii_ - q2i_) T'' " ^ (aipi + fci_L, {a - ai)pi + fc2_L, -api - fcij_ - /c2±) 



We have omitted terms ~ I3p2 in the arguments of d^,, since they do not contribute to d^^^d^^ , see Eq. ( 10 ). Introducing 

2 j,2 ^ 

the variable u = ai/a and taking residues at /3i — and /JJ^ = we obtain 



9 rani rbn' l' 



d^zd^z' / d^qid^q2d^kid^k2 eK9i+92,^)x-i(fei+fe2,a)i-»(9i-fci,^)-ife-fc2,^') 



(gi± + g2.L)A(fci±-Hfc2±)A' dK(""Pi + gi-L)^V(^«Pi + ki±) d^^juapi + q2±)d\,{uapi + k2i_) 
(91 + '72)i(fci + ^2)! ql^u + ql^u kl^u + kjj^u 

X V^^^iuapi + qi±,uapi + q2±, -api - qij_ - q2±) T'' " ^ {uapi + ki±,uapi + k2±,-api - ki± - k2±) 



where we have imposed a cutoff a < in accordance with Eq. ( 17l 
Using formulas 



df_,^{uapi + g_L)d^^ {uapi + k_i_) = (^g^^ 



2 

-(91 + 92)ir^^A(a'"Pi + qij_,aupi + q2j_, -api - {qi + q2)±){g'^ - — P29i)(3± - 7£:H^'2'72) 
= iq!± ~ <&i}9l + ^91(91 + 92)^ - \{qi + q2y4 
we can represent the contribution of diagram in Fig. Isk in the form 
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du uu I d^qid^q2d^kid^k 
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(23) 



2 2 1 r 2 2 

(9? - 92)'^ij — qu{qi + q2)j + -{qi + 92)192^ (fc? - ^2)'^^ — fcii(fci + ^2)^ + -(fci + ^2)1^2^ 

u u JL u u 

Throughout the paper we use Greek letters for indices ji — 0,1,2,3 (with 5^'^=( 1,-1, -1,-1)) and Latin letters for 
transverse indices i — 1,2. 

The diagram shown in Fig. K5p is obtained by the substitution e^**^'^i+'"'^'^^'' _g-«(fei+fe2,2;)i (^^^ic different sign 
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FIG. 4: "Cut vertex" diagrams. 



comes from the replacement [~oopi,0]y by [0, —oopi]^). We get 



(24) 



du uu I qid^ q2d'^ kid^ k2 



J{qi+q2,x)±-i(qi-ki,z)-i{q2-k2,z') |'g-i(fci+/c2,2;)j p-i(fcl+fc2.' 

{qi + q-zVih + hViqfu + g|u)(fcfS + kju) 

2 2 1 r 2 2 
(9i - ^D^ii qii{qi + Q2)j + -{qi + 92)4(72^ {kf - kl)5rj + k2)j + -{kx + k2)ik2j 

U U J L It U 
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B. "Cut vertex" diagrams 



Next, consider the "cut vertex" diagram in Fig. 4r. The analog of Eq. (20 1 has the form: 



dt du dv{A'i{tpi + x±)Ai{upi+yi_)A'',{vpi + y^_)) 

Jo J — 00 J u 

= ig^sr""" j daid(3ida2dP2 j (f 'gi* "gs^'fei^'^a T^'^aipi + <7i±, a2Pi + g2±, -(«! + "2)^1 - (gi + (72)^) 
ai0(ai)0(a2)e»(«^+«--)- ((gi + 92)^ + 2(/3i + /32)P2) ^ (^i + ^-^^^P^) ^ (^2 



2(fc2,g2)j 



(ai + 02) (/3i + /32 - «e) 



[(ai + 02) (/3i + /32)s - (gi + 92)1 + ie] 

^-i(ki+k2,y)±Jjmbjjnc^-i{qi-ki,z)±-i(q2-k2,z')± 

kj{kfa2 + k^ai){aif3is — qfj^ + ie){a2l32S — q^^ + ie) 



(25) 



Going to variables a = ai + a2, u ^ ai/a and taking residues at /?i + = and ^ we get 
g^J dt I du I dv{Al{tpi + XA^)A'l{upi+yj_)A'',{vpi + yj_)) 

-,4 



oo J u 



fj^ r<^ dry r r pi{qi+q2.x)±-i{ki+k2,y)±-i(qi-ki,z)±-i(q2-k2,z')± 

2^2^ io «'io J J {qi + q2)l{qlu + qju)kl{klu + kju) 

2{ki,qi)±_ 



X (fci + ^^^P2) (fc2 + ^^2,_g2j± ^ \ ^ g2)±Ar^'^\attpi + 91^, a^lpi + 92±, -api - {qi + 92)±) (26) 
V aus / II V aus / v 
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and therefore 



(Tr{C/,C/t})Fig|4], (27) 



/O " 



(<?!+ 92)^(9?" + g|u) ukl{klu + klu)^ 



where we have used the formula 

aus / aus / v 

r 2 2 1 

= kiik2j {ql-ql)5ij qu{qi + q2)j + -{qi + q2)iq2] (28) 



u u 



following from Eq. (|22|). 

The contribution of the diagram shown in Fig. Wb differs from Eq. (27) by the substituion e"**^''!^'^^''')^ 
g-i(fci+/c2,a;)j_ changing the order of , matrices. (Similarly to the case of the Fig. diagram, this prescription 
follows from the replacement [— oopi,0]y by [0, cx^pij^; but now we consider the second term of the expansion in the 
gauge field). We get 

(Tr{f7,?7^t})Fig|4].+b (29) 

= -^^2 f - tduUu[d^k,ct-k2t\,d\2!d^zd-z'UT'U-r {ql~ql)k,-lqu{q, + q2), + l{q^ + q2U2, 

271-2 7o a Jo J J {qi + q2y{qfu + qiu) 

ukf{klu + k2u) ^ ^ 

There is another type of diagrams with two gluon-shockwave intersections shown in Fig. [5] 

/■oo />0 /"O 

dt du dv {Al{tpi + Xi_)Al{upi + xj_)Al{vpi + yi_)) (30) 

J J — 00 J — 00 

= 2g\s j daida2d;Pid;p2j d^zd^z'V^'^Uir j '^^^^^J^^ j d^qid;^q2 e{ai)e{a2) 

X r^'T^^^^aipi + qi,a2Pi + 92, -(«! + a2)pi - (91 + ^2)^) e^l-Ji.— -)^+*(92,x-^')^-<('ci:X--)^-»(fe2,y-.')x 
[{qi + q2)± + 2(/3i + P2)P2]X fciM + ^^^5f^P2^ A:2. + ^^^^P2. 1 



(ai + a2)(/3i + P2- ie)[{c(i + a2){Pi + ^2)8 - {qi + 172)1 + ie] aif3is - + ie a2/32S - (?|_|_ + ie 
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Taking residues ai Pi + P2 — and at /3i = and going to variables a ^ ai + a2 and u — ai/a we get 



Jo J — 00 J — 00 

^i{qi,x — z)±+i{q2,x—z')±—i(ki,x-'z)±~i{k2.y~z')j 



9^ ramn /""^ 

2^2^ Jo a^Jo 
V sau / V 



2(92,^2)^^2 



9 ramn 



01- Jo 



kfkliqi + q2)^{uqj + uq^) 
{qi + q2)^Tfj,„x{aupi + qi,aupi + q2,-api - {qi + q2)j_) 

i{qi ,x — z)±+i(q2,x—z')±—i(ki,x—z)±—i{k2.y—z')± 

kfkliqi + q2)^{uql + uq^) 



{ql± - ql±W' - -qliqi + q2y + ^{qi + 92)'?^ 
u u 



kuk2j 



and therefore 



(31) 



(Tr{C/,;7t})Fig[5] 



' da 











du / d^zd^z' / d'^qid^q2 / d^kid'^k2 e-*(«i"'=i'^)^-»fe-'^2^^')-LC/™"' 



J(qi+q2,x)±-i(ki,x)±-i(k2,y}± 

{qi+q2)^{qlu + q^u) 



{ql - 91)% - -91.(91 + 92)j + 4(91 + 92).92, %§^Tr{i"f/,f"'r'f/t} 
u u i k^k2 



where again we have used formula ( 28 1 . 

The sum of the contributions ( 24 ) , ( 29 ) and (l32|) can be represented as follows 



(Tr{C/a;;7^})Fig[3]fFig|4]fFig|5] = (Tr{;7a;C/j{})Fig|6] I+III+V+VII+IX 



A_J^^J^ du uu j d^zd^z' jd\id\2 ji^kid^k2 e-(9i-fci^-)^-%2-fc2,.')^[/Wc/?"' 
i(qi+q2,x)±_ r . ^ 2 2 



X /"^"'"Trlt" 



(91 + 92)^(9?" + 9iw) 



(9i - 92)% - y9i-*(9i + 92)i + ;^(9i + 92).92i 



Ux 



, , ( p-i(ki+k2,x)^ p-i(ki+k2,y)_L\ y 

b fbm 11 \^ ) 



(fci + k2)'^{kju+ k^u) 



{kf - kl)5ij fcii(fci + k2)j + -(fci + k2)ik 



'2i 



/ ^n' -j-m' ^-i{ki+k2,x)± g-iiki,x) ±-i{k2.y) ± , , ^m' ^n' ^-i{ki+k2,y) ± 

4:iki k2 I t^^ + 

* "'I ?2fc^(fc^-u + fc|w) uuklk^ uki{klu + k2u) 



-1} 



(32) 



(33) 



If we add contribution of the diagrams with the gluon on the right side of the shock wave attached to the Wilson line 

p*(9i+92,y)i) obtain 



at the point y instead of x (which differs from Eq. (33 1 by the substitution e*'^''i+'^'^)^ —e" 



(Tr{f7a; J7j{ }) Fig|6] 1+11+. . .+x 
S'^'^Jo a Jo 

(pi{qi+Q2,x) 



(34) 



du uu d^zd^z' / d^qid\2 / d^kid^k2 e-'(<li->'i,-)^-<<l2-k2,z')^jjrara'^nn 



ai{qi+<!2,y)A 



(91 + 92)^(91'" + 9|m) 



(9i - qD^i] ~ ~'?i»(9i + 92)j + 3(91 + 92).92j 



Ux 



( „—i{ki+k2,x)j p-'iiki+k2,y)±\ 

b rbm n K'^ J 



2 

{kj - kl)5ij - -ku{ki + k2)j + j{ki + /c2)ifc2j 



2 

-( 

u 



Aik\ik 



(fci + k2Y{klu + klu) 

^n' -j-m' ^-i{ki+k2,x)±_ g-i(ki,x)±~i(k2,y)± , , ^m' -j-n' ^-i{ki+k2,y) 



liK2j I 



ukf{k1u + fc^w) 



uuklk. 



ukl{kfu + fc|u) 



The result ( 34 1 can be obtained from the self-energy contribution ( 24 1 by the replacement of the term corresponding 



to the emission of the two gluons via the three-gluon vertex 



b rbm 



(e 



— i{^k\+k2,x) A 



-i{ki + k2,y)A 



(fci + fc2)2(fc2u + fc|u) 



(fci - k2)6ij - -kii{ki + k2)j + 3(^1 + k2)ik2j 
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with similar contribution containing the "effective vertex" 



S™"(A:i,fc2;a;,y) = 

( p-i{ki + k2,x)j p-i{ki + k2,y)±\ 

b fbm 11 \^ ^ l_ 

(fci + k2Y{klu + klu) 

^n' -^m' ^-i(ki+k2,x)± g~i{ki ,x) ±-i{k2 ,y) ± 



2 2 

[kl - kl)5^, kii{ki + fca),- + -(fci + fc2)ifc2, 

u u 



Aikiik2j 



^ f7yi' -j-n' -i(ki + k2,v) ± 
l"' _|_ 

ukf{kfu + klu) uuk-fk^ ukf{klu+k2u) 



(35) 



It can be demonstrated that the sum of the contributions of the diagrams shown in Fig. [6] I,..., IV, XI,..., XVI can 
be obtained from the self-energy contribution ( 24 1 by replacing the gluon vertex 



(pi{gi+q2,x}± _ pi{qi+q2^y)±\ 
. 'i^ J_ 

{qi + q2yiqlu + qju) 



with similar "effective vertex" 



a j:amn 



i{qi+q2-x)A pi(<?i+<?2,j/)i 



Cj-m-i-n ^i(qx+q2,x) ±_ gi(qi,x) ±+i(q2,y) 
uql{ql 



2 2 

[ql - ql)5tj - -qiiiqi + q^)] + ^(91 + 92)»92j 



2 2 
[ql - 92)% — qiMi + q^)] + -(qi + q2)tq2] 



(36) 



uuqlql 



infm^i{qi+q2,y)± 

uqf{qfu + q'^u) 



(37) 



Note that (37 1 is equal to S'^'^^{qi,q2;x,y). Let us consider now the box diagrams topology shown in fig. [6]XVII- 
XXXIV. The calculation of these diagrams is similar to the above calculation of "cut self-energy" and "cut vertex" 
diagrams so we present here only the final result 



( Tr { J7:e C/y } ) F ig |6] X VI I + . . . +XXXI V 



(38) 



du uu / d^zd^z' / d'^qid^q2 / d^kid:^k2 e-*(«i-'=i'^)-L-^fe-'==.^')^t/f ™'c/^r' 



X Tr< 



27r2 7o a 7o 

/ ^m^n^i(qi+q2,x)± ^n^m ^i(qi+q2 ,x) ± ^i(qi,x) ±+i{q2,y) i_ 



^^'^^M uql{qlu + qju) 



uq^{qlu + qju) 



uuqlql 



-ee 



a'i{qi,y)±+i{q2,x)_i 
uuqlql 



fn^m ^i{qi+q2 ,y) ± frn-fri ^i{qi+q2,y) ± 



uql{qlu + qlu) uq^qfu + q^u) 



kiik2j 



^n' -^m' ^—i{ki+k2.x)± ^m' ^n' ^—i{ki+k2 -x) ± ^—i{ki,x)^—i{k2,y) 

uk'l{k\u + k^u) ^ ukKkfu + klu) uukfkl 

^m'^ji'g-i(fci+fc2,y)i iri' -j-m' ^~i(ki+k2,y)±_ 

^ uk\{k\u+klu) ^ ukl{k\u+klu) 



-i{k2,x)±-i{ki,y)± 



This expression agrees with the sum of "box topology" diagrams in Ref. [TT] . 

Now we observe that each three-gluon vertex diagram is equal to its own cross diagram (the same cannot be said 
for box diagrams). Thus we may redefine the "effective vertex" (351 in the following way 



5™" ik,,k2;x,y) = 



b £hra n 



i2kiik2j 



{ki + k2)'^{klu + klu) 
ukl{klu + fc|u) 



2 2 

{kl - kl)5ij - -kii{ki + k2)j + 3(^1 + fe)ifc2i 



-i(ki,x)±-~i{k2.y)± 



-i{k2,x)±~i{ki,y)± 



ukKkfu + k^u) uukfkl 

^m' -j-n' ^-i{ki+k2,y)± fri -^m' ^-i(ki+k2,y) ± 



ukKkfu + k^u) ukKkfu + k^u) j 

which corresponds to writing each contribution of the three-gluon vertex diagrams as a sum of two equal terms. 



(39) 




FIG. 6: Diagrams with two cuts. 



12 



A similar expression can be written for the "effective vertex" (37) and therefore the sum of all diagrams with two 
gluon-shockwave intersections can be written as 



871-270 a Jo 

X C/r'"'c/;,"'Tr{[5t"-('Zi,'72;x,y)C/J [5"'"' (fci, ^2; y)f/t] } 
Separating the contributions of different color structures one obtains 



(40) 



(41) 



87r2 







d^zd^z' / d^qid;^q2 / d^kid^k2 t/f ™ U^, 



t 



(^^i{qi,x)+i(q2,x') „ 2; ^ |-(^2 „ ql)uu5ij - 2uqu{qi + 92)^ + 2u((7i + 92)492^ '7h92j . -'7i»92j 



{qfu + qju) 



(91+92) = 



9? 



X 



9i 



(/c^u + k^u) 



2i 



_ jbni'W klik2j ^^-i(kuX)-i(k2,Y') _ ^ ^ 

This result agrees with Ref. [11] 
Performing the Fourier transformation 



{ki + k2)^ 

y) - z{t™',r'}%^(e-^('=-^) - e'^('=-^))(e-^('=-^') - e-^('=-^'))}c/, 

1 2 



we get 



d'q,d\2 e'('--^)+*(«--^) f , , 

qfiqtu + q^u) 



XiiX2j 



[ d^qid\2 e'(9i-=^i)+*(9^^^2) " - igH(gl + g2)j- + |(gl + g2)»g2j- 

y ^ ^ (gl +g2)2(g?U+g|u) 

_ -(X^ - X^)(5tj + Ijxi - X2)tX2^ + Ixiijxi - X2)j 

47r2(a;i — X2)'^{uxl + uxl^ 



(42) 



— (Tr{C/.C/t})p,,,6, 
X Tr{iV 



r X. 



rabc 



XiX' XiX', 



(z~z')^ ^ uX^ 



[tX'2 



+ uX'^ 

where we introduced the notations 

— / v2 



du UU I d^zd^z' Ufuf 



{x^y) 



l^a' ja'6'c' 



MuX2y'2 



X^Y' 

- (x^y) 



uuXW^ 



^ uu VX2 y2;^;^,2 Y'^Ji 

uu Vx2 r2;^;^,2 Y'^Ji y ^ ' 



= {X'' -X")5,j + -{z-z'),X'^ + -X,{z-z')j 
Y,, = {Y^ -Y'^)6,, + \z-z%Y;+'^Y,{z-z'), 



(44) 



C. Subtraction of the (iO)^ contribution 



It is easy to see that result (43 1 for the sum of diagrams in Fig. [^diverges as w ^ and u ^ 1. If we put a 
lower cutoff a > ct' on the a integrals we would get a contribution ~ In ^ coming from the region 0:2 ^ ai > a' (or 
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ai ^ a2 > u' ) which corresponds to the the square of the leading-order BK kernel rather than to the NLO kernel. To 
get the NLO kernel we need to subtract this {LO)^ contribution. Indeed, the operator form of the evolution equation 
for the color dipole up to the next-to-leading order looks like 



d_ 
d-q 



(45) 



where rj = \na. Our goal is to find i^NLO by considering the l.h.s. of this equation in the external shock-wave 
background so 



d 



Shockwave — Shockwave 



Shockwave 



The subtraction (46) leads to the ^ prescription for the terms divergent as ^ (and similarly ^ 



contribution divergent as u — s- 1 ) . Here we define 



in the usual way 



du f{u) 



11 _ r,, /w-/(o). 

u 



u 







du f{u) 



11 _ fiu)-m 

u 



u 



+ 



(46) 
for the 

(47) 







To illustrate this prescription, consider the divergent terms in Eq. (43) proportional to {X,Y){Y' ,z — z') or 
{X',Y'){Y,z-z') 



drj 



aj. du 
4^ 



^fa'b'c'-j-a' 
L*y2 + uY'^ ^ Y^ 



d'zd^z'ufu^f (^'^)(^''^-^') 



(z- z')2y'2 



1 



iX^ + uX' 



{t''',i^'})c/t} + .r'^'-'Tr{( 



+ uY'^ Y^ 



72 + ^O^nk.^'^'c^J} 



-^^^-^^—^^-^riir'^e + i{t'x})u,r'''^'e'ul + er^''^u^{r''''''e' ~ i{t^' ,t^'})ul} 



,.2 /•! 



du 

47r4 Jq u 



d^zd^z'ufU7'^-^^f^^^ 
{z~z'YY^ 



1 

liX2 + uX' 



uf ^ ^ t" 

{uY^ + uY''^) ' Y 



72T + 7^0^^^'})c^.n + r'^'^'Tr{(-f 



uYiYU 



I -pabc 



lY^ + uY"' Y 



-+'{t'',n)u.yui} 



/"^^Tr{-r[/,( 

+ ^,2, J . Tr{(-^''-t°+^{^^t-})£/.r^'-'^°vt-^v°'^^.(^'''^'^°'+^o'^^^'})t/,^} 

X' (uY'' + uY' ) 



(48) 



Note that the second term is equal to the first one after the replacement m <-*■ m, z ^ z' and b c, b' c' . 

It is convenient to return back to the notation ai and a2 — (t — ai (after the value of a is set equal to a) 

d 



' '"'da2 f,,,2jrrbb'rrcAX,Y){Y',Z^z') 



d^zd'z' C/f V:? 

' {z-z'YY'^ 



(49) 



r^-Trlt^C/J — 



ah' c 4-a' 



aiX^ + a2X'^ 



r2 + a2Y'^ Y^ 



iiY^ + a2Y' 



2ia 



X^iaiY"^ + a2Y'^) 



The corresponding term in i^LoTr{J7j,C/^} is (see Eq. Hi) 



4a, 



d^z 



(50) 



(a;-z)2(y-z) 

The relevant term in the "matrix element" (-ftrLoTr{J7j;C/|}) in the external shock-wave background comes from Ux, Ul 
taken in the leading order in a, (so that ^ U^, Ul U^) and 11^ (8) taken in the first order in a. 



{U^ <E> ill) 



a, da 



Jo a.2 



{z - z'Y{y - z')"" 



(51) 
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or vice versa: Ux ^ Ux, Uz ^ Uz and 

a 



TI" Jo 



[z - z'Yiy - z'Y 



(52) 



Here we have used the leading-order equations for Wilson lines with arbitrary color indices [UlIS]- Substituting eqs. 
( [51) and ([52]) in Eq. (|50| we obtain 

2al rda2 /■^2^^2^, (^'^)(^''^-^')^„r,,,aVc',c,fc,,,6 



{Ki^oTriUxUh) 



TI" Jo OL2 



X^Y^Y'^{z- z') 



-j:^i:v{{if'''^'''' tH'^Uxt^' Ul - if^H^'Uxt^'f'uDufu'^'} 

(53) 



From Eq. ( 46 1 we get 



(i^NLoTr{C/,;7]}) 



a1 rda 



27r4 7o a2 

pah c -la 



, ^ -fa b c fi 

f-'^-Tr{t''Ux(— - ^ 

^ \aiY^ + a2Y'^ Y^ 



2al rda2 i ,2 



{z^z'YY'' 
{t^',i^'})c/t} + Tr{( 
2i(T 



(54) 



^ + ^{^^n)f/.t'''^/t} 



X2(air2 + Q,2y. 



TT Jo 0:2 J X^Y^Y' [z - z'Y " " 



2^7o T 



d^zd^z'vfVT (^'^)(^''^-^') 



(z-z')^>"' 



r''7°'^'"'Tr{i'^t/^t°'[/|} 



2u 



y2 



1 



1 



iX2 + uX'^ x^ 



Y2 



2i 
X^ 



1 1 

uX^ + uX'^ ^ X^ 
1 



{uX^ + uX'^){uY^ + uY'^) X^Y^ 
r'''''Tr{{t',nUxt^'ul} 



lY^ + uY''^ Y^ 



Tr{tH''Uxr'''''''t'''ul ~ tT'"=Uxt''' f^' U^} 



which corresponds to the 



prescription (47 1 (the same prescription was used in Ref. [TP). Note that the "plus" 



prescription (47 1 is a consequence of the "rigid" cutoff \a\ < a (17 1; with the "smooth" cutoff ^ we would get 
different results - see Appendix B. 



D. Assembling the result for 1— >3 dipoles transition 

There are four color structures in the r.h.s. of Eq. (|43|). Three of them reduce to 



r^^r'''''''ufu^'Tv{eUxe'ul} = iTr{C/,C/t}Tr{[/,t/],}Tr{C/,,[/t} - i Tr t/^t/.C/],} + [z ^ z') 

^f-bc^fu-/Tr{t-Ux{t''\t''}Ul) = - iTr{[/,;7l}Tr{C/,C/],}Tr{[/,,C/t} + iTr{C/,C/t[/,,C/t [/,[/],} - (z ^ z') 

ir'^'^-'ufu^Tiiit", t-}Uxe'ul) = iTr{[/,C/t}Tr{C/,C/],}Tr{C/,,C/t} + ^-Tv{UxUlUz'UlUzUl} - [z ^ z') (55) 

We will not need the explicit form of the fourth color structure U^'^ U^f' Tr({i", t''}Ux{t'^ , jC^jJ) since it is multiplied 
by pure LO^ integral /|^ = /^ + /17 and does not contribute to the NLO kernel. 
Performing integration over u using the prescription ( 47 1 after some algebra we get 

»i 



du 







1 



iA2 + ?2A'^ V( 



X 



1-3 



uA2 



X,Y^ 



{z-z'Y 

VA2r'2 -x'2y2 



A2r'2 + A'V2-4A2(z-z')2 A2r'2 
- 4 + 2 ^ ^ In 



A2y'2 - x'2y2 



1 



1 



A2r'2 Y^x' 



{x-yY 

{z - z'Y 



1 



1 



A2r'2 x'2y2 



In 



A2y^ 
A'^r2 



(56) 



15 



and 



du 



X,, 



In 



X^Y'^ 



'^'nx^J 



-{x^y) 



1 



uuX^Y'^ 
1 











VX2 " 


Y2j 







X2Y'^ X'^Y^ 



In 



X^Y'^ 
X'^Y^ 



so the two-cut contribution ( 43 1 reduces to 

„2 



drj 



X2y'2 _ x'^Y^ 
r (x-y)^ 



167r4 



1 



d'^zd'^z' 

1 



X2Y'^ Y^X'^ 



I (z-z')4 

(x-j;)2r 1 



X'^Y'^ + X'^Y^ 



Aix-yfiz-z') 



M2 



(Z-Z')2 



1 



^2y/2 J)j^/2y2 



In 



X^Y'\ 
X'^Y^i 



[Tr{U,Ut}Tr{U,Ul,}Tr{U,,Ul} - Tr{U,UtU,,UlU,Ul}] 



{x-yf ( 1 



In 



^2y/2 



Z — Z' )2Vx2y'2 y2;^/2;/ x'2y2 
This result agrees with the 1^3 dipoles kernel calculated in Ref. ^TT]. 



Tv{U^Ul}Tv{U,Ul,}Tv{U,,Ul} 



(57) 



(58) 



E. Subtraction of the UV part 



The integral in the r.h.s. of Eq. (58 1 diverges as z ^ z' . It is convenient to separate the divergent term by 
subtracting and adding the contribution at z = z' : 

Tr{C/,(7l}Tr{{/,;7l,}Tr{;7,.C/t} - Tv{U^UlU,>UlU,Ul}] (59) 
= [Tr{t/,C/l}Tr{C/,C/],}Tr{C/,,[/t} _ Tr{;7,C/]C/,-C/tC/,C/],} - [z' ^ z)] + [Tr{[/,C/t}Tr{C/,C/t} - Tr{C/,C/t}] 

For the last line in the r.h.s. of Eq. ( [58| the subtraction is redundant since 

{x - yf ix-yf f 1 



X^Y'^ Y^X 



In 



^2y,2 



X'-^Y2 







Jf2y/2^,2y2 (^z - Z' 

The easiest way to prove this is to set y = and make an inversion x ^ \/x so the integral ( 60 1 reduces to 

(i — z)2 



,2 ~l {x Z, X z') 

d z —■ — TT In ■ 



{x~z'Y{z-z'Y (x-z'y 



Thus, we obtain 
d 



drj 



(Tr{C/,(7t})p.^l^ 

(x " y)^ 
X2y'2 - x'2y2 



167r4 



d^zd^z' 



( " (z - zM4 + 



1 1 

X2y'2 ^ y2j^.2 



4 

(z ^0'^ 

(^ - yr r 1 

(Z-Z')2 



j^2y/2 ^ ^,^Y'2 - 4(.T - y)2(z - z'f 

(z-z')4[x2y'2 -x'2y2] 







^/2^2 



■In 



X2y'^ 



X'2y2 



J 



+ 



(a; - y)'' 



X2y'2 X'2y2J 

X [Tr{C/,C/t}Tr{[/,[/],}Tr{C/,,C/t} - Tr{C/,C/t[/,,C/t _ (/ ^ z)] 
" C ( ^2^^ + ^) } 1- fSTr{f/.C/l}Tr{L/.C/t}Tr{f/.C/t} 



a: 



167r4 



X^Y'^X'^Y^ (z-z')2Vx2y 



d2z[Tr{i7,!7l}Tr{C/,t/n - Tr{(7,C/n] / d^z' 



(60) 



(61) 



(a^ - y)"^ 

X2y'2 - x'2y2 



1 1 

X2y'2 ^ y2x'2 



{z~z'f 



1 



(z-zO 

1 



X2y'2 x'2y2 



} 



X2y'2 + X'^y2 - 4(a; - y)2(z - z') 



^'2^2 



M2 



(z-z')4[x2y'2 -x'2y2] 



■In 



x2y^ 
X'2y2 



(62) 
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The first term is now finite while the second term contains the UV divergent contribution which reflects the usual UV 
divergency of the one-loop diagrams. To find the second term we use the dimensional regularization in the transverse 
space and set d± = 2 — e. Because the Fourier transforms (42) are more complicated a,t d± ^ 2 it is convenient to 



return back to Eq. (40 1 and calculate the subtracted term in the the momentum representation. The calculation is 



performed in Appendix A and here we only quote the final result (164) 

4 r^X^Y'^ + X'^Y^ -4{x-y)^{z- z')^ 



drj 



(Tr{[/,f7j})Fig[6] 



167r4 



1 



a zd z 



V {z-z'Y 
^ vf r 1 

{z-z'Y 



{z - z'Y[X^Y'^ - X'^Y^] 



1 



X^Y' 



J X'^Y^J 



X [Tr{[/,C/t}Tr{C/,C/],}Tr{C/,,t/t} - Tr{(7,C/lC/,.C/tC/,C/],} - (z' ^ z)] 



1 



;i^2y/2^,2y2 [z - z'yKx'^Y''^ Y^X' 

67_^ 
9 3 



2 /2 

} f7^^Tr{t/,C/l}Tr{C/,C/],}Tr{C/,,[/t} 



{x-vYiU^^ X^Y^ 
3 (a; - yY 



X2Y2 



[Tr{C/,C/t}Tr{C/,C/t} - _Tr{C/,C/t}] 



where /i is the normalization scale in the MS scheme. 



(63) 



IV. DIAGRAMS WITH ONE GLUON-SHOCKWAVE INTERSECTION 

A. "Running coupling" diagrams 

The relevant diagrams are shown in Fig. [t] (plus permutations). Let us start from the sum of diagrams Fig. [7]a 
and b. It has the form: 



du / dv {A'^iupi + x±_)A,{vpi + y±))Fig^+h 



(64) 



"da f ^ , fd(3df3' 



{k^ + ieY {k-k'Y + ie 



r^'^^((a - a')Pi + {k~ fc')±, a'pi + k'^,-api - k^) 



{k + 2pp2)y 



xT''''''^'{{a~a')pi + {k-k')j_,a'pi + k'j^,-api-kj_) - 2 



{k + 2/3p2).(g,, + "^^2,0 g^'dlik') - d^^-ik') 



(fc2 + ieY 



k'^ + it 



where the first term in the square brackets comes from Fig. [7^ and the second from Fig. [7]3. We use the principal- value 
prescription for the 1/a' terms in dp,u(k') in loop integrals. 

To regularize the UV divergence we change the dimension of the transverse space to 2-e. After some algebra one 
obtains 



FIG. 7: "Running coupling" diagrams. 



du / dv {A'^iupi + xj_)Al{vpi + 2/_L))Fig[7}i+b 



g^Nc^J."^^ J d^-^kd^-^k'd^-^q J d^-^z UfJ^ da J d(3dp' - 



^i{q.X)^-i[k.Y)i_ 



X - / da'dfi' 



X <i - e[{a - 2a')(3s + k'l - {k ~ k')l][ „ (fc, g)± + {2k' - k, g)^] + 2^ ^(fc,g)i/3s 



a 



{a - 2a') 



a 

a + a' 

a — a' 



(fc, q)^{2k' - fc, fc)_L + (a - 2a')(3s{2k' - fc, q)x_ + (g, fc)±(fc - 2k', k - fc')_L + (g, fc')±(fc, 2fc' - fc)_ 



a - 2a' 



[q, k)^{k, fc - fc')± + 2fc' - fc)i(fc, fc - fc')± + (g, fc)±(fcl - ki) + {q,k- k')^{k, 2k - fc') 



+ (a - 2a'){q, k - k')^f3s + {q,k~ k')^{k'\ - (fc - k')\) - {q, k')^{k', fc - fc')± + {q,k- k')^{k ~ fc')i 



a' — 2a 
a' 



a — 2a' 



a 



(g, fc)i(fc, fc')i + (g, 2fc' - fc)i(fc, fc')± - (g, fc')±(A:, k + k')^ + (g, k)^{k', 2k - k')^ 



+ (a - 2a'){q, k')^l3s + (g, fc')±(fc, 2fc' - fc)^ - (g, k')^k'l + (g, fc - fc')±(fc', fc - fc')± 
a + a' a' — 2a 



a — a a 



(fc, fc')±(g, fc - fc')± + (g, k')±{k, fc - fc')± + [(fc - fc')i + fc ](g, fc)± + 



4a(g, fc)j 



a — a a 

— + 7 

a a — a 



(3s — 4Q!(g, fc)j 



/2 



r fc' 



a ~ a' 



(fc - fc') 



'^2 
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where we have omitted the contribution 



1 



a'{a'(3's - /c'l + ie) 



1 



1 



k'^ + it ■ '{k',P2) 



(66) 



Taking residues at /3 = and = and changing to variable u = ^ we obtain 

o M 

du / dv {Aliupi + Xi^)A\{ypi + y_L))Fig[7]i+b 

^ — oo 

' ' ' a Jo kS^[k'\ + {k - k')^u] 



(67) 



8vr2 J J J, 

X I - 2£(2fc' - k, k){q, k' - ku) + (1 - 2u)(fc, q){2k' -k,k) + (q, k){k - 2k', k - k') + [q, k'){k, 2k' ~ k) 
l + u 



u 

2-u 



(1 - 2u){q, k){k, k^k') + 2{q, k)k'^ - 2(fc, fc')±(g, k') 
(1 - 2u)((7, fc)(fc, k') + 2(g, fc)(/s - fc', fc') - 2(g, fc')(fc - fc', A:) 



(1 + m)(2-w) 



[(fc, k'){q, k-k') + {q, k'){k, k - fc')] + {q, k) (fc - k'f + fc'^ - 4 



fc'' .(fc-fc') 



l\2 



Using the "plus" - prescription (47 1 to subtract the (LO)^ contribution we get 
'^-(Tr{C/,C/t})Fig.[7],+b 



drj 



(68) 



2a27V^^2e /^2-e^[Tr{f/,f/t}Tr{(7,C/t} - --Tr{U.,Ul}] d^-^kd^-^k'd 



Nr. 



1-2-e^ 



k^q"^ 

(q.2k-k') (q,k + k')\^ [k-k'f f\ (2-e)(q,ku~k')(k,k~2k') + 2(q,k)P 

— '- T r, In T 1- du ^ 

fc'2 {k-k'yj fc'2 7o fc2[fc'2e+ (fc- fc')2u] 

Next we calculate diagram shown in Fig. [Tj:. 



POO /"O 

<?M dW du / dv {Al{tpi + x)i^(upi + y)A',{vpi + y))Fig[7l= 

•/ •/ — oo -/ ti 

= -~2g*f'"'J dadpda' dp' dp" d^-^qd^~'kd^-^k' e^'''^''''^^^'^''^y-''^^ 



(69) 



g(a)L/°' (fcl + 2/3V2)M 

(/^(a/Ss — fc2 + ie) a'(3's — fc'2 + ie 



{q± + ^(g, fc)±P2);, ((fc - fc')^ + 2p"p2), r^"^(aVi + fcj, (q - gQpi + (fc - k')^,-api - kj_) 
^ a'{a-a') {a - a')P"s - {k - k')l + ie {(3' - ie){P" + f3' - ie){P - (3' - f3" - ie) 

There are 2 regions of integration over a's: a > \a'\ and a < \a'\. Taking relevant residues, we obtain 



i(q,X)~i{k.Y) i-l 

d'^-^kd'^'^k'd'^~^q I d^z±_ Uf ^ du , ., 

' " k'^q^k^ io U'^fi+(fc-fc')2u 



.Jl^fabl^2e I da I ^2-ei,^2^e,J^2~e„ / ^2 ^ , ttcI ^ 



lk'^u-\-(k - k'] 



i^[k"+k'] 



•J2 



+ — (g,2fc-fc')-2(g,fc)(fc',fc-fc') 



1 



(fc - k'Yu + fc^i 



-[(g, fc)fc'' + (g, fc')fc'] - (fc', 2k - k'){q, fc)] I (70) 
u J J 



where we have introduced the variable u = |a'|/a as usual. After integration over u with help of Eq. (47 1 this reduces 
to 



-l-^f-bi^2e £ ^ j ^^-^kd^-'k'd^-'q I d^-'z^ uf 



^^^^[(g,fc)fc'' + (g,fc')fc'] + (g,fc) 



In 



^i^[fc'Vfc^] + (,,fc'-2fc) 
2(g,fc)(fc',fc-fcO | 



In 



(fc - fc')= 



(fc - k'f r^^^ 2{q, k){k', fc - fc') j e^(9^^)-*(fe''*') 
P Jo " fc'^u+(fc-fc')2uJ fc'\2p 



(71) 
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and therefore 



dr] 



^{q.x)-^{k,Y) ,r (q,fc') (g,fc')fc2 (g,A:'-2fc)i (fc - fc')^ 



(g,fc) , , fafc-fcQfc^]^^ P ^ r'^^ {q,k ){k',k~k' 



k'^ 

'■'A \ f72) 

' ' k''^{k-k'yi^" k'^ "Jo^"' k'^[k'\+{k-k'yu]i 

Next we calculate the sum of diagrams in Fig. [7]d,e, and f. The contribution of the diagram shown in Fig. |7jl is 

{Tr{U,Ul})p,,,^ (73) 



aa'iP + 13' - ie){P - ie)(a'/3's - k']_ + ie){a(3s -kl + i€)ql 



" Jo 



where we took residues ai f3 — —, f3' — —+ and introduced the variable u ~ ° , . It should be noted that the 
cutoff a < (7 in the r.h.s. of this equation translates into dada' 6{a — a — a') while our cutoff ( 17 1 corresponds to 
J^dada' 9{a — a)9{a — a'). Fortunately, the difference 

f°° dada' ,, ^, 1 dv ^ fc'^ + fc^w 

/ 7, \ei(j ~ a)e{(j - a') - e{(j - a- a')] = — In (74) 

does not contain In a and hence does not contribute to the NLO kernel. Similarly, one can impose the cutoff ai+a2 < o 
instead of the cutoff ai, Q!2 < cr in other diagrams whenever convenient. 

Before calculating the diagrams in Fig. [7^ and Fig. [Tf it is convenient to make the replacement 

du dv dt A''{u)A''{v)A''{t) -> - / du dvdt A''{u)A''{v)A''{t) (75) 



2 



which can be performed since the color indices b and c in ...t^f^... are contracted. For the diagram in Fig. [t]; we get 
(Tr{J7,C/t})Fig[7^ = 2gV'"cfTr{i"f/,t^f/t} (76) 

. j,a,a',p,p'j,^--k,^-^k',^-^qjd^-^.Uf ^^^'^"^^-^'^^"^^ ( , P' \Oia)ik,q), 



'{13 - te)2{k'^ + te) \p + P'-ie /3 - - leJ a{k^ + ie)ql 

= 4gV^'cj.Tr{rC/:,i''C/]} J dad:a'd:(3dp' J d^-'kd^-'k'd^-'q J d^-'z ;7afcg^(«,x)^-^('c,i')x 

^ 0{a){k,q)^ 

^ aa'{(3 ~ ie){a'f3's ~ k''\ + ie)(/3 + /?' - ie)(/3 - (3' - ie) {a/Ss -kl + ie)ql 

where cp = ^yn~- Taking residues at /3' — —(3, (3 — ^ aX a' > Q and /?' = /3, /3 = ^ at a' < we obtain 
(Tr{C/,[/t})pig.[7], 

drv r i{q,X)±-i{k,Y)± f f ( h rA , 



The diagram in Fig. |7f yields 

(Tr{J7,C/t})Fig.[3 (78) 

= oV^^ / dad^d^-^k I da'dp'd^-^k' ^ 

J J a'{f3-ie){f3' -ie){a'p's-k'l + ie) 



27r2 Jq a Jq uuj J J k'^k'^q^. ^ 
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Adding Eqs. ( |73| , ( |77| , ( |78| and integrating over u using Eq. (|47| we get 

d 



(Tr{C/,C/t})Fig.[7,j+,+f 



(79) 



Note that the diagram in Fig. [7f does not contribute to the NLO kernel. 

The contribution of the last "running coupling" diagram shown in Fig. [7^ has the form 



(80) 



oo pO 



= - ^Ti{t''Uj''t''t'^Ul} / du / dv{A^{upi + x±)A'i{vpi + yj_)) / dtdw{A\tpi + x±)A''{wpi + .t_l)) 



where we have again replaced j'^^dt p_^dw A^{t)A'^{w) by ^J^^dtdw A''{t)A'^{w). Using the Eq. (jllj) we get 



{Tr{uM})Fi^.m 



(81) 



i9_ 

-,4 



PI PI 



da da' 



\x) 



Pj^ Pj^ Pj^ Jo a a' 

which is obviously a (LO)^ term which does not contribute to the NLO kernel. 
Combining expressions ( 68 1 , ( 72 1 , and ( 79 1 we get 



dri 



(Tr{C/,;7t})Fig.[7] 



(82) 



^,,(g,2fc-A:')-^^^(a,fc + fc') In 
{q,k') , {q,k')e , {q,k'~2k)^^^{k~kr 







\ (2 - e)(q, ku - k')(k, k - 2k') + 2(q, k)k^ 
du 2 



k'u+{k-k'Y 



u 



{k-k'Y k'\k-k'y 



k'' 



{q,k) ^ {q,k) ^ {q,k-k')P^^^k^ f^^^ {q, k){k' , k - k') ^ 2{q,k) ^^k^ 



(fc-fc')2 k'\k-k')^ 



k' 



k'^[k'^u + {k- k')^u] 



k'^ k'^ 



1 r i(q.,X)^l(k,Y) { l2 

2alN,^l^' I d^~'z[Tr{U,Ul}Tr{UMl} - — Tr{C/,C/t}] j d^-'kct^-'k' d^-' q ^-^ <^ In ^ 



3(g,fcOfc^ ~4(g, k){k,k') 
^ k''^{k-k'Y 

Using the integral over k' 







\ {2-e){q,ku-k'){k,k-2k') + 2{q,k)k'^ f\ 2(q,k)(k',k ~ k') 
du ^ — du — ^ — ^ 



k'^[k'u+{k - k')^u] 



/c'^[fc' u + (fc-fc')^"] 



^2-efc/ j 3(g,fcOfc2-4(g,fc)(fc,fcO ^^ fc^ 



I k''^{k-k'Y 

2{q,k){k',k-k') ] 
10 k''^[k''^u+{k-k'yu] I 

one reduces the r.h.s. of Eq. ( [82| to 



-(Tr{C/,f7t})Fig.[7] 



drj 

X fd^-^kd^-^ 



\ (2 - £)(g, fcu - /fc')(^, - 2fc') + 2(g, fc)fc2 
du 2 



fc'-' 7o k'^[k'u+{k-k'Yu 
(g,fc)/ r(£/2)r^(l-|) rll ^^2^ 



47r 



1 (/c2)e/2 r(2-e) 



2^ ' 



q e 



1 d2--z[Tr{;7,C/]}Tr{[/,C/t} - -IxrlC/.L/t}] 

^5} 



(g,fc) f r(£/2)r2(i-|) rii y 

^ ^ ' 3 6 



Pg2 l(fc2)£/2 r(2-e) 



(83) 
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Next we subtract the counterterm 



22 alNc 



Tie 



(84) 



corresponding to the poles 1/e in the loop diagrams in Fig. |7](we use the MS scheme). We obtain 



:iq.x)^i{k,Y) (g, fe) I'll 



k'^q^ 



3 k^ 



67 



3 i 



(85) 



The complete set of running-coupling diagrams is presented in Fig. |8] 

The contribution of diagrams in Fig. [s] VII-XII differs from Eq. ( |85| by the replacement e'^'''^-' e^iiX) ^nd 
sign change. There is also a symmetric set of diagrams XIII-XXIV obtained by the reflection of diagrams I-XII with 
respect to x* axis. The result is obtained by the substitution e"*^'^'"^) Q-^{k,Y) therefore the contribution of 
all diagrams in Fig. ([8| takes the form 



X Jd^kd^k'd\ [e»(9.^) -e^(9^^)][e 



d^z[Tr{U,Ul}Tr{U,Ul} 



1 



^ k^q^ 13 k^ 9 3 



3 i 



(86) 



The remaining diagrams XXV-XXVIII contribute only to the (LO)^. We have shown this for diagram XXVII (Fig. 
[7^), see Eq. (81 1. The diagram in Fig. |8] XXV is obtained from the equation (81 1 by the replacement x ^ 
y, and the diagrams XXVI and XXVIII by the replacements {x\p-U^\y){x\:^\x) by (a;|^[/^|x)(y|^|y) and 

{y\^U^\y){x\^\x), respectively. Thus, the diagrams XXV-XXVII do not contribute to the NLO kernel. 

There is another set of diagrams obtained by the reflection of diagrams shown in Fig. with respect to the 
shock-wave line. Their contribution is obtained from Eq. ( |86| b y the replacement q ^ k in the logarithm so the final 
result for the sum of all "running coupling" diagrams of Fig. 18] type has the form 



dr] 



(Tr{C/2;[/j{})Fig.|8] total 



2n 



d^z[Tr{U,Ul}Tr{U,Ul} - ^Tr{[/,C/t}] 



(87) 



87r3 



X J d^kd^k'd^q [e 
''d'z[Tr{U.,Ul}Tr{U.Ul} - ^Tr{C/.[/t}] 



-i{k,X) 



^ fc2q2 1 3 q2fc2 



2^1 X^Y^ 

^ In :r- 



134 
~9~ 



2^ 
3 



1 11 


■ 1 


1 ■ 


J + y 







134 27r2 
~9~ ^ ~3~. 

^2. 



In 



r2 



B. Diagrams for 1^2 dipoles transition 

There is one more class of diagrams with one gluon-shockwave intersection shown in Fig. [9] These diagrams are 
UV-convergent so we do not need to change the dimension of the transverse space to 2 — e. First we calculate diagrams 
shown in Fig. |9k,b. 



(Tr{C/,C/t})Fig.|9j 



+b 



ATr{eUj''t''t''ul} / d^-^z / daidl3ida2d(32 / d^^^kid^~''k2 e 



-i(fei ,y)±-iik2,y)j 



Sbcijad 

(/32-«e)(/3i+/32-ie) 



e{ai) / (j;2^^iiq.x^z)^+i(kuz)^ 



{I3i-ie){l3i + l32-ie) 

gi{k2,x)± 



ai{ai(3is - kfj_ + ie)q\ Q:2(a2/32S - fc|_|_ + ie) 



9 



"da 



du / d^^'kd'^^'k' e-*('='«)^-^(''-''2')^ 



Sbdljac 



^abjjdc 



k'^u{k'^u + k''^u) k'^u{k'^u + k'^u) 



^2^^i(q,x-z)j_+i{k,z)j_+Hk' ,x)j_ (g' 



^TT{t'^U,thH''ul} [ d^z H— [ d^kd^k'd\ ^ 
71-2 " J Jo a J g2 



sbcTTad ^bdrjac 1,2 
(q,X)~i(k,Y)^ + i(k' ,x-y)^ a -0 K_ 

k^k'"^ k'^ 
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(XXVI) .' (XXVII) (XXVIII) 



FIG. 8; The full set of "running coupling" diagrams. 



and therefore 

^(Tr{;7,C/t})Fig.[9}.+b (89) 
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FIG. 9: 1^2 dipoles transition diagrams. 



The contribution of diagrams shown in Fig. is obtained from Eq. (89) by the replacement a; <-> y in the left part 
of the graph and the sign change so that Q-'^Hy-^)+tik .^-y) ^ _^-ik(x^z)-i{k ,x-y) ^ ^j-^^ smn of the diagrams in Fig. 
|9^-d takes the form 



^(Tr{J7,;7t})Fig.|9,,_d = d'z[Tr{U,Ul}Tr{UMl} - ^Tr{C/,[/t}] 



i 



i{q,x-z) ( -i(k,y-z)±-i{k' ,x-y)± 



a; <-^- y ) In 



(90) 



Next relevant diagram is shown in Fig. [9fe 



OO J —OC' 



g" I du I dv I dt {A'^{upi+x)A,{vpi+x)Al{tpi + y))Fig^ 



= 2g^J d:aidPida2dP2ct^kiS'^ki 



ejai + a2)f'''{ki^ + 2/?iP2M)(fci + W2P2.) 
aiQ!2(/3i - ie){ai(3is - kl^ + ie){(32 - je)(a2/32S - kjj^ + ie) 



X / d^z U. 



ad 



oi{q-ki,x-z)±~ik2{y~z)± 



2(fcl+fc2,g)j 
(q!1+C(2)s 



P2 



g2[(ai + a2){Pi + f32)s ~ {ki + /ca)! + ie] 



'-T^''\k^,k2,-k^~k2) 



(91) 



There are three regions of integration over a's: ai,a2 > 0, ai > ~a2 > and a2 > —a\ > 0. Going to the variables 
a = ai + a2, u — a2/a in the first region, a — ai, u — —a2/a in the second and a — a2, u — — ai/a in the third, 
we obtain 
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/■oo (.0 /-O 

„3 / J„, / J„. / J-L I AO./ I ^\ Ab / 



oo — oo 



g-^y du / dv I dt {Al{upi + + x)Al{tpi + y))Fig[9}= 



(92) 



'da 



27r2 Jo a Jo 



du / d'^kii'^ki f^'^uf 



J(q-ki,x-z)±-i(k2,y-z)j 



uki, kt + 2 



''Ifj. \ '^21/ 

2 



2(fci + fc2,g)± \ r^'^^(aMpi + kj.aupi + fc^, -api - fcj^ - k^) 
as P'^^J uukfj_k2j_{ki + ^2)^^ 

(A:i + /c2)i \/ 2(fci + fc2,g)_L \r^''^(api + fcjL^-Mapi + fc^,-au-fcjL-A;J-) 
= P2u 9a H = P2\ 



ukf^iki + k2)]_[u{ki + /C2)i + ukj 



2±J 



(ki + 2^-'^l±J^p^ k-L 
\ aus J u V 



2(fci + fc2,9)± \ Tf^^^l-uapi + ki-^api + kk^-au-k}- -ki) 
= P2 



A uk'^^{ki + k2)\[u{ki + k2)\ + ukl 



i±J 



Using the formula 
2A 



+ — P2) (ki + — P2) (q+ — P2) T^"'^{aipi + ki^,a2Pi + k2±,-{ai+a2)pi-{ki + k2)i_) = 

\ S / fi\ S J u\ S /A 

-C(ai - a2)(fci, /C2)_L - A{ai + 2a2)(g, k2)± + B{2ai + a2){q, ki)± - [{q, ki)±{k2, ki + k2)± - {ki ^ fe) 



2C 



we get 



(Tr{J7,?7t})Fig.|9j 



87r2 



cf^A:i(f^fc2^ / d^z 



_^ (1 + u){ki + fc2)^(g, fci) - (1 + u){q, ki + fc2)(fcl, ^2) - fcl)(fc2, fcl + k2) - {q, fc2)(fcl, fcl + fc2)] 

'Ufc2(fci + k2Y{klu + (fci + fc2)2u) 

_^ -(1 + U)(fcl + fc2)^(g, k2) + (1 + fci + fc2)(fcl, fc2) - M[(g, fcl)(fc2, fci + fc2) - (g, k2){ki,ki + fc2)] 

ufc|(fci + k2)^{kju + (fci + fcz)^^) 



(93) 



Performing the integration over u (with prescription (47l) we obtain 
d 



dr] 



(Tr{C/.[/t})pig.|9^ 



9*K 
87r2 



^i(q-ki,X)-i(k2,Y)jjc 



f (fci + fc2)^ [(g, fc2) In ^^^^ - (g, fci) In ^^^^^^l - (g, fci + fc2)(fc? + fc2') In % ] 



y"d2z [Tr{[/,C/]}Tr{t/,C/y^}-Tr{[/,J7]}] j d^kd^k'*^ ^l{q-k\x-z)-^{k-k' ,y-z)jjc 



87r2 

(g,fc-fc') 



I- 

Ifc' 



In 



k''^{k-k'f fc'2 k''^{k-k'Y {k-k'Y k^k 

where we made the change of variables fci — > fc' and fc2 — > fc — fc'. 
The sum of diagrams shown in Fig. |9k-e can be represented as 



(g,fc)^^(fc-fcr 



(g,fc) ^^ (fc-fcQ^ - 



fc2(fc-fc')2 



_d 



(Tr{C/,?7t})Fig.[9j,_, = 2a2iV,y" d2z[Tr{;7,C/t}Tr{C/,C/t} - -lTr{C/,C/t}] ^ d^kd^k'd^ 

{q,k-k') , ((z,fc) (g,fc)- 



X e 



i{q,x-z) ( -i{k,y~z)±-i{k' ,x-y)± 



X -f-^ y 



k'^{k-k'Y k'^ik-k'Y k^k' 



In- 



(94) 



(95) 



Note that the expressions (90 1 and ( pM| are IR divergent as fc' ^ but their sum (95 1 is IR stable. Once again, the 
contribution of the diagrams in Fig.T9t-k are obtained by replacement e^'^^^~^^ —^My-^) go the contribution of 
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diagrams of Fig. [9]a-k has the form 



drj 



{q,k-k') 



{q,k) {q,k) 



-k'^{k-k'Y k^{k-k'y k^k'^ 
Performing the Fourier transformation with the help of the formula 



In- 



(96) 



a:kdk 



^ g-t{k,y)-i(k' ,x-y) _ 



kj k 



Va;2 



{k~k'Y fc2 k^{k-k'Y 

2 „,\2 



In 



^'^ln (^7^ In^^ 



» ( {x,y) 

sttA 



(here n = yjx-^ip — (x, y)"^) one obtains 
d 







)-\ 




1 IK 








)-\ 




1 IK 





In \ 



In 1 



U — ^ "2 U — ^ 2 



In u 



In u 



W— M ^^T9 U — ^- 



1 In — In (^~^'^) 

2 {x-y,y)-iK 



- In ^ In ^ 



(x-yY 



{x,y) - 



C?77 
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(97) 
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Note that the two last terms in the r.h.s. of Eq. (97) do not contribute. 

The contribution of the diagram obtained by reflection of Fig. [9]with respect to the shock wave differs from Eq. ( 95 1 
by replacement q k which doubles result ( 98 1 . The final expression for the contribution of all "dipole recombination 
diagrams" of Fig. [9] type has the form 
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V. ASSEMBLING THE NLO KERNEL 



Adding results ( 58 1 , ( 87 1 and ( 99 1 one obtains the contribution of the dia, 
with the shock wave in the form: 



^ (Tr{ }> Fig[6lf Fig[8lf Fig|9] 
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{x - yf ^ {x - yf / I 1 \ -1 . 

72 j^/2i 

mere are also diagrams witnout giuon-snocKwavc mterscction iiKe tne grapn sn^ . ^. r--~c- 

to the parent dipole 'Yt:{UxU^} and their contribution can be found from Eq. (100 1 using the requirement that the 
r.h.s. of the evolution equation must vanish at a; = y (since UxU^=l). It is easy to see that the replacement 



^ + ^ : „ + ' „ ) \ In — ^Tr{UxUl}Tr{U,Ul,}Tr{U,jUl} 

X^Y' X' y2 (z — z 

There are also diagrams without gluon-shockwavc intersection like the graph sho wn i n Fig. 10 They are proportional 
to the na,rent dinole TtUJ„U1\ and their contribution can be found from En. (100 
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FIG. 10: Typical diagrams without the gluon-shockwave intersection. 



Tt{U^UI}Tt{U,UI} - ^Tt{U^UI} by Tt{U^UI}Tt{U,UI} - iVcTr{ C/^C/t} fuifiUg the above requirement so one 
obtains the final gluon contribution to the NLO kernel in the form 
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Promoting Wilson lines in the r.h.s of this equation to operators and adding the quark contribution from Ref. [3] 
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(102) 



we obtain the full NLO kernel cited in Eq. ([5| 
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VI. COMPARISON TO NLO BFKL 



A. Linearized forward kernel 



In this section we compare our kernel to the forward NLO BFKL rcsults[13 . The hnearized equation ([5| has the 
form 
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(103) 



For the case of forward scattering (JJ{x,y)) — {U{x — y)) and the hnearized equation (103) reduces to 
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Using the integral J13 from hep-ph/9704267 [TS] we get 
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drj 
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B. Comparison of eigenvalues 

To compare the eigenvalues of the Eq. ( |106[ ) with NLO BFKL we expand U{x, 0) in eigenfunctions 

e»<^(xiM^)^ (Zi(n,7)) , (107) 



(Z^(xi,0))= ^ I '^'"^ .»<Ar^2„2^7 



n— — oc> 2 



— — zoo 



27ri 

compute the evolution of (t/(n, 7)) from Eq. ( |106| and compare it to the calculation based on the NLO BFKL results 
from [13 HH]. (For the quark part of the NLO BK kernel the agreement with NLO BFKL was proved in Ref. [Tf\). 
The relevant integrals have the form 

„2 



^ 'd^z [2(z2/x2)^e™^-l] ^ 



27r 



(x — z)2z2 



= x("-,7) 



i/<i-|2(.V.^)V.*- 11(^4) 



1\ (x-z)2 47x(7) 
— I In 5 = X (", 7) - X ('^^ 7) 



^2 _ "1 
' 4 



d^z (zVx^)'^ 



m0 



(x — z)2z2 



In 



(x — z)2 z2 

^ ^2 2 



X^ x^ 



In^ = ;^x"('^,7) + x'(",7)x("-,7) 



where x(",7) = 2(/'(l) - "0(7+ f ) - V(l - 7 + f ), and 



2x2z2 f 2 2 X^ - ^2 



16x2z2 
4 , 0^4 o^2^2 



Z^ 



(108) 



(109) 



00 



dt „ ^ 1n-^ + ^ 



x2+t222 

2 



77 



iV3y2(3-27)(l + 27)"'"; (1 -27)sin2^7 



I — 

i fi 



cos 7r7 



F(n,7) 



1 



(x2 — z2) r x2 x2z2(x — z) 

In — In 



2^ 7 (^Va;')^ 'e"*<^ _ ^^2(2; + 2)2 L"- ^2 (a;2 + z^Y 



2L.2(-^j-2L^-- 



(x2 - z2)2 



(x — z)2(x + z)2 



du v^z"^ 
In 



where TSl 



$(n,7) 



(x — Zu)2 x2 



<i>(n,7) - <i>(n, 1 - 7) 



1 + t 



{y^-2V/(— )-Li2(t)-Li2(-t) 



(0(. + 1) - m + Ml + ^) + E ^) - E 7fc^[i - (-1)'^]} 



fc=i 



(110) 



(111) 



29 



The convenient way to calculate the integrals over angle (j) is to represent cos ncj) as T„ (cos (j)) and use formulas for 
the integration of Chebyshev polynomials from Ref. [18] . 

Using integrals (108) - (110) one easily obtains the evolution equation for U{n,^) in the form 



{U{n,-i)) = \ 
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where x'(n,7) = ^x(",7) etc. 

Next we calculate the same thing using NLO BFKL results |T31[IH]- The impact factor $a((z) for the color dipole 
U{x, y) is proportional to as(g)(e''^ — e''^)(e"*'^'^ — e^*"?^) so one obtains the cross section of the scattering of color 
dipole in the form 
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where G^{q, q') is the partial wave of the forward reggeized gluon scattering amplitude satisfying the equation 

ioG^{q,q')^5^^\q-q') + j d'pK{q,p)G^{p,<^) (114) 
and <^B{q') is the target impact factor. The kernel K{q,p) is symmetric with respect to q <-!■ p and the eigenvalues 
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The corresponding expression for {U{n,j)) takes the form 
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where 9 is the angle between q and x axis. Using Eq. (114) we obtain 
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The integration over q can be performed using 
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(recall that K{q,p) = /ir(p, q) and as{p) = ag ~ In ^ with our accuracy). The result is 
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FIG. 11: Gluon loop inside the shock wave . 



where the angle <p corresponds to p. Since w ~ as we can neglect terms ^ uj m the argument of i5 and expand 
x(n,7— ^) ~ x("-j7) ~ f x'('^i7)- Using again Eq. (1141 in the leading order we can replace extra w by ^Ncxi^,^) 
and obtain 



d as nn r(-7+f) [d'^pd^q 
s— (6/(71,7)) — — -^r^ COS — I — —e 



ds 



2 r(i + 7 + f)y p2 ^,2 - vv 

X $5(9') / "^^°°^(^)''G.(p,g')^^c[xK7) - ^X'K7) + ^['5^7) - 2x(n,7)x'(n,7)]] (120) 



Finally, expanding a'j,{p) ~ as{p){as - ^ In ^)as{fx) we obtain 
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which can be rewritten as an evolution equation 
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X" (n, 7) - 2x(n, 7)x'(«, 7) + 6C(3) + F{n, 7) - 2$(n, 7) - 2$(n, 1 - 7)] } (^(^, 7)) 



(122) 



This eigenvalue coincides with Eq. (1121 up to the extra term 2^(3) . It would correspond to the additional contribution 



to the r.h.s. of eq. (5 1 in the form of "1^2" C(3)Tr?73;£/j which contradicts the requirement -^jUxU^ — a± x = y. A 
possible reason for tne disagreement is the connection between the matrix element of the color dipole with a rigid 
cutoff a < a and the cutoff by energy s in Eq. (1131. It is worth noting that the coefficient 6C(3) in Eq. (1221 agrees 
with the j ^ 1 asymptotics of the three-loop anomalous dimensions of leading-twist gluon operators [19] . 

It should be emphasized that the coincidence of terms with the nontrivial 7 dependence proves that there is no 
additional 0{as) correction to the vertex of the gluon - shock wave interaction coming from the small loop inside the 
shock wave, see Fig. [ll](In other words, all the effects coming from the small loop in the shock wave are absorbed in 
the renormalization of coupling constant in the definition of the U operator (|6|). In the case of quark loop, we proved 
that by the comparison of our results for Tr {UxU^} in the shock- wave background with explicit light-cone calculation 
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FIG. 12: Renormalon bubble chain of quark loops. 



of the behavior of Tr{UxU^} as x ^ y [Qj. For the gluon loop, we can use the NLO BFKL results as an independent 
calculation. Let us repeat the arguments of Ref. [S] for this case. The characteristic transverse scale inside the shock 



wave is small (see the discussion in Ref. [9j ) and therefore the contribution of the diagram in Fig. 11 reduces to the 
contribution of some operator local in the transverse space. This would bring the additional terms with the nontrivial 
z dependence to the kernel which translates into the nontrivial additional 7-dependent term in the eigenvalues. Such 
terms do not exist and therefore the gluon interaction with the shock wave does not get an extra 0(as) correction. 



VII. ARGUMENT OF THE COUPLING CONSTANT IN THE BK EQUATION 



In this section we briefly summarize the results of the renormalon-based analysis of the argument of the coupling 
constant carried in Refs. [Oj [lOj 

To get an argument of coupling constant we can trace the quark part of the /3-function (proportional to n/). In 

2 

the leading log approximation In ^ ^ 1 , <C 1 the quark part of the /3-function comes from the bubble chain 
of quark loops in the shock- wave background. We can either have no intersection of quark loop with the shock wave 
(see Fig. 12 1) or we may have one of the loops in the shock-wave background (see Fig. 12 3). 
The sum of these diagrams yields 



dr] 
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aic/"^(/) 



(123) 



where we have left only the /3-function part of the quark loop. Replacing the quark part of the /3-function —f^n-f In ^ 

2 

by the total contribution j^bln^ we get 



iip,y)A 



,,^^^as{{p-lf) 



(124) 



J p^ - ■ ■ ■-■ [p-iy 

In principle, one should also include the "renormalon dressing" of the double-log and conformal terms in Eq. (|5|. 
We think, however, that they form a separate contribution which has nothing to do with the argument of the BK 
equation. 

To go to the coordinate space, we expand the coupling constants in Eq. ( 1241 in powers of as = as{iJ?), i.e. return 
back to Eq. (1231 with fjJ^/ — > —bf^- Unfortunately, the Fourier transformation to the coordinate space can be 
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performed explicitly only for a couple of first terms of the expansion as{p ) ~ — ^ Inp /fi + Inp //i ) . In 
the first order we get the running-coupling part of the NLO BK equation ([5| 
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The result of the Fourier transformation up to the second order has the form |S] llOj 
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(126) 



We extrapolate the In + In terms in the above equation as follows: 
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(128) 



where dots stand for the remaining conformal terms and In^ term. (Here we promoted Wilson lines in the r.h.s. to 
operators) . 

When the sizes of the dipoles are very different the kernel of the above equation reduces to 
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F - y| < F - z\, \y - z\ 
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(129) 



In the earlier paper[2j the Eq. (1271 was interpreted as an indication that the argument of the coupling constant is 
the size of the parent dipole x — y. We are grateful to G. Salam for pointing out that the proper interpretation is the 
size of the smallest dipole as follows from Eq. ( 129 1. 

It is instructive to compare our result to the paper [lOj where the NLO BK equation is rewritten in terms of three 
effective coupling constants. The authors of Ref. [TD] extrapolate Eq. (1261 in a different way 
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1 /y2^j. 1 (v2. 2{x-z,y-z) a,(X2)a,(y2) 
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where is some scale interpolating between X^ and (the explicit form can be found in Ref. fW). Theoretically, 
until the Fourier transformations in all orders in Inp^ /^^ are performed, both of these interpretations are models 
of the high-order behavior of running coupling constant. The convenience of these models can be checked by the 
numerical estimates of the size of the neglected term(s) in comparison to terms taken into account by the model, see 
the discussion in Refs. |22j 



VIII. CONCLUSIONS AND OUTLOOK 



We have calculated the NLO kernel for the evolution of the color dipole. It consists of three parts: the running- 
coupling part proportional to /3- function (see diagrams shown in Fig. [s]), the conformal part describing 1 — )■ 3 dipoles 
transition (diagrams in Fig. |6| and the non-conformal term coming from the diagrams in Fig. (|9|. The result agrees 
with the forward NLO BFKL kernel [13 up to a term proportional a^C{3) times the original dipole. We think that the 
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difference could be due to different definitions of the cutoff in the longitudinal momenta (see the discussion in previous 
Section) . It would be instructive to get the j — > 1 asymptotics of the anomalous dimensions of gluon operators directly 
from Eq. ^ , without a Fourier transformation of our result to the momentum space and comparing to NLO BFKL 
as it is done in Sect. |VTj The study is in progress. 

There is a recent paper [20j where the dipole form of the non-forward NLO BFKL kernel is calculated using the 
non- forward NLO BFKL kernel|21]. The kernel obtained in [20 is different from our result (and not conformally 
invariant). We think that at least part of the difference is coming from the fact that the evolution kernel ^ should 
be compared to the non -sym metric "evolution" NLO BFKL kernel K'^™^{q,p) rather that to the symmetric kernel 
K{q,p) defined by Eq. (113). The kernel K°™^ corresponds to the Green function G^^ defined by Eq. ( 113| with 



different lower cutoff for the longitudinal integration 

1 f(Pa(fa' f'^~^^°° duj / s \^ ~ 

{U{xM - ^2j^^-sme^^^-l){e-^''^-l)^BW)j^_^^ G.(,,g') (131) 

The G^{q,q') satisfies the equation (1141 with the kernel i^"™' 



u:G^{q,q') = S^^Hq~q') + J d^pK'=^°\q,p)GUp, q') (132) 
and the relation between K°^°^{q,p) and K{q,p) has the form (cf. Ref. [13 ) 

K'^^°\q,p) ^ K{q,p) - \ j £q'K{q,q')\n^K{q' ,p) (133) 



2 



It is easy to see that the structure ( 131 1 repeats itself after differentiation with respect to s so it can be rewritten as 
an evolution equation for U{x) (whereas the derivative of the original formula ( |113[ ) does not have the structure of 
the evolution equation due to an extra j^)- In terms of eigenvalues, the modified kernel (1331 lead to the shifts of 



the type x('^:7) ~^ x(^:7 ~ which we saw in Sect. VI 3 



It should be emphasized that the conformally invariant NLO kernel describes the evolution of the light-like Wilson 



lines with the "rigid" cutoff in the longitudinal momenta (17). On the contrary, for dipoles with the non- light-like 
slope the sum of the diagrams in Fig. [6] is not conformally invariant (see Appendix). The reason is that a general 
Wilson line is a non-local operator which is not conformally invariant to begin with - for example, the non- light- like 
Wilson line turns into a circle under the inversion — > x^/x^. With the light-like Wilson lines, the situation is 
different. Formally, a Wilson line 

[oopi + Xi_,—'X)pi + x±\ = Pexp <ig dx^ A+{x^ , x±) > (134) 

J —oo 

is invariant under the inversion x^ x^ jx^ (with respect to the point with zero (-) component). Indeed, (x+,a;j^)^ = 
—x\ so after the inversion xj^ x^lx\ and x'^ —> x^ lx\ and therefore 

f 2:+ a;+ ~| 

[oopi + — cx)pi + Fexp <ig d— A+{ — ,x±)> = [cx)pi + — oopi + a;j_] (135) 

^ J— 00 

Thus, it is not surprising that the bulk of our NLO kernel for the light- like dipoles is conformally invariant in the 
transverse space. The part proportional to the /3-function is not conformally invariant and should not be, but there is 

another term ~ In ^^"^ja In \y~^Jj2 which is not invariant. The reason for that is probably the cutoff jaj < cr which can 

be expressed as a cutoff in longitudinal coordinate x^, and therefore under the inversion x^ x'^ jx\ the cutoff can 
pick up some logs of transverse separations. It is worth noting that conformal and non-conformal terms come from 
graphs with different topology: the conformal terms come from 1— >3 dipoles diagrams in Fig. (pi which describe the 
dipole creation while the non-conformal double-log term comes from thel— >2 dipole transitions (see Fig. |9]) which can 
be regarded as a combination of dipole creation and dipole recombination. It is possible that in the effective action 
language, symmetric with respect to the projectile and the target |23j . the evolution kernel is conformally invariant. 
We hope to study this problem in a separate publication. 
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IX. APPENDIX A: UV PART OF THE ONE-TO-THREE DIPOLES KERNEL 



As we mentioned above, it is convenient to separate the UV-divrgent and UV-finite parts of the Eq. (43 ) by writing 
down u^"^'UIT' = (^i/nim' jj^m' _ f^mm'jjmm') _,_ jjmm'jjnm'^ rj.^^ Contribution of the first part leads to Eq. ([58 1 
while the second UV-divergent term have the same color structure as the leading-order BK equation. After replacing 
Ijrnm jjnn j-^y ijmm jjnn ^ integrating over u with the prescription (47 1 and changing variables to k2 = q2 = k' , 



p — qi + q2, I — qi ~ ki (so that qi — p — k' , ki ~ p — I — k' and ki + k2 = p — I) the Eq. (43 ) turns into 
dt^-'pa:^-' F,{p, l)+ j i^pctH F2{p, I)] (e*(P^^) - e^(f'^))(e-^(P-'-^) - g-^^P"''^)) 



where 



X 



+ (e- 



-i{p-l,X) _ ^-i{p-l 



J{p-k\X)+i{k\Y) _ i{p~k',Y)+i{k', 



.f\2 



(k',p~k')(p^k'Y -2(p-k',p-l-k'){k',p-l-k') , ip-l^k') 
X In 7\ 



{p-lY{p-k'Yk'\p-l - k'Y 
+ j d^pd^W^k'{{e'^P^''^ - e^ip^Y)^(^e-KP-^-k',x)-^(k'.Y) _ ^-^ip-i-k' ,Y)-^ik' ,x)^ 

(k',p- I - k')(p - k'f - 2(p - k',p - I - k')(k',p- k') , (p - k'f 

In ■ 



p'^{p-k'Yk'^{p-l-k'f 



(136) 



1-f J ^ {p-k'f + {p-k' -If {p-k'f k'^ + {p-k'f {p-k'f 

pHp - ir 1 {p - k'Y -{p- k' - lY ""{p-k'- If ^ (p - k'Y - k'^ 



{p- k' -l)^ + k'^ ^^{p~ k' -If ] , 2{p,p-l)lf {p-k',p-k' -I) 

[p-k' - ly - k'^ ^ k^ J 



p^{p-iY \{p-k'Y{p-k' -ly 



1 \ , (p-k'f (p-k' -If 

In 1 



( 



{p-k', p-k' -I) {p-k',p-k'-l) ^ ln{p-k'f/{p-k'-lf 



2 {p,p-l-k') p2 



(p - k'f 



{p-k' - If 



V {p - k'f -{p-k' -If \ p^{p-l- k'fk''- 



■In 



2 {p-l,p-k') ^^ {p-lf 



{p-lf {p-k'fk'' 



(137) 



and 
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Pi , {p-k%\ \n{p-k'f/{p-l-k') 



p2 (p - A;')V A;''[(p - k'f - {p - k' - If] 

/2 7/^2 



(138) 



1/fc" 



In 



{p-lf {p-k' -lyJ {p-k'Y -{p-k' -ly {p-l-k') 



2{l,p-k' -l){p, k')/k 



{p-k'f (p-k\p-l-k'){p,k') ^ {p-k'f 
■ in -, — r-TT ^ m 



{p - k' - iY[{p - k'Y ~ {p - k' - ly] {p~k'~iY k'^{^p-k'Y{^p-k' -ly k'^ 



+ 



+ 



+ 



{i,k')/k 



/2 



{p-k'Y {p-k'-l,k') ^^ {p-kr I {p-l,k')/k'^ ^ Jp-k'-iy 



{p-k'Y-{p-k'-if'{p-i-k'Y 2k'^{p-i-k'Y (^p-i-k'Y-k'^ 

2 {p-k'f 2ln{p-l-k')yk'^ {p,p-l-k') {p-k'f 

{p-k'y-{p-k' -If {p-i-k'Y {p-i-k'Y-k''' ^ {p-i-k'Yk''' p2 

-2{l,p- k'){p - I, k')/k'^ {p-k'f _ {p-k',p-l-k'){p-l,k') {p -k'- If 



ip-iy 



{p-k'Y[{p-k'Y-{p-k' -ly] [p-k'-iy k'\p-k'f{p-k' -If 



k'^ 



il,k')/k 



jp-k'f ip,k')/k" {p-k'f ip-k',k') {p-l-k'f 



+ 



{p-k'f-{p-k' -If {p-l-k'f {p-k'f-k'^ k'^ 2k'^{p-k'f k'^ 

2 {p-k'f 2\n{p-k'f/k'^ ^ {p-l,p-k') ^^{p-l-k'f 



{p-k'f-{p-k' -If {p-l-k'f (p-k'f -k'^ ip-k'fk'^ {p-lf 
We need to perform the integration over k'. Let us start with the UV-divergent term ^ Fi. Using the integrals 



(139) 



k''^{l-k'f k' 



r(d-i) ;4-d 



47r / (t'^k 



In -TT In 



, {k',k'^l) ^^^{p-k'f l^^p\^{p-l) 
k'^{l-k'f 

, (p.k',p-k'-i) ^ {p-kriP-k'-if ^ _^p^]^{p-i)' 



2 P P 



f{p-k'-lf 



k'^ 



P P 



A.U'k' iP'P-i-k') ,^ P' 



{p-l- k'fk'^ k'' 



r(|)r(f-i)r(2-f) 

= {p,p-i) 



(140) 



one obtains 



Flip, I) 



I r 2ip,p- I) r^(l - I) f _ 4 + i^U '^^'^ " In ^ 

t7r[p2(p_|)2 r(2_£) ^ / >\ ^ p^{p-lf\2. 



47r 



- in g in - in- + -g) - - ^^^^ ,7 1^^/^^ + Oje) 



P P 



P 



'2 ' 3 / 3(p - If 3p2 



(141) 



Let us at first consider the UV-divergent contribution 

v2 



d 
da 



(p-i.jc) _ ^-i(p-i.r)^ {p,P-l) [ ^-(1 2) 



p'^ip-lf I r(2-e) 



^r(./2)(-4 + y) + 



(142) 
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To this contribution we should add the counterterm corresponding to quark and gluon loops lying inside the shock 
wave. The rigorous calculation of the counterterm was performed in Ref. [9J and the result is 



da 



{Tr{uM}) 



a? 2 



CT 



,N, 



1. 



- / d'-'z±{^Tr{U,.Ul}Tr{U,Ul} - -Tr{C/.C/t}) 



J p^{p-iY 



(143) 



where we need the gluon part of h {— ^Nc). After subtraction of the counterterm ( 143 1 the UV-divergent contribution 
( 142 1 reduces to 



da 



{Tr{U^Ul})vy- 



CT 



p^{p — lY 



11 , P 67 

— In — r 



so one obtains the regularized Fi in the form 



FT^{p,l) 



1 

47r 



67 



p2(p-;)2V3 9 



p2 



4n' 



lnpVZ2 ln(p-Z)2//2 



i{p~lY 



3p2 



It is convenient to calculate first the Fourier transform with e*^^''"^) '■'^p i-Y) ^ Using the integrals 



d^pdH e*(P'^)+'(''^) (f;^~/) ln4 
d'^v dH e»(P'^)+*(''^) ^P^P~^^ In ^ In ~ 



1 ^ 



47r2 X2r2 

2 



A2 



1 (x,Y) , ^2, r2 

In In 

47r2 X2y2 A2 A2 



d'p d^l e 



p^ip^lf 



1 {X.Y) 
47r2 ^2^2 



In" 



X2 

r2 



we get 



1 



167r3 X2r2 



Y2 

2 ^ 



2X^ 



1^ A^+^''' y2 



27r2 



67 

33 



487r3 



1 y2 1 x2 

:^i^^A2 + 72i^A2 



Hereafter we use the notation A = X — Y = x — y. 

Next we calculate the F2 contribution. We need the following Fourier integrals: 



/ 



(p-l), {p-l-k% 



Pi 



{p-k% ^ Injp-k'f/jp-l-k'Y 
{p - k'Y ) fc'2[(p - fc')2 -{p-k' - If 



1/fc 



,2 



In 



{p-k') 



J\2 



{p-iy {p~k' -ly) {p-k'Y -{p-k' -i)'^ {p-i-k'Y 



{X,Y) 
167r3X2y2 



ln2 



X2 

y2 



2^ 



(144) 



(145) 



(146) 



(147) 



(148) 



d^p dH e'(pA)+.(/,i') 
' Uk' 



{l,p-l-k'){p, k') {p - k'f 

p2 J k'^p -k' - if{{p - k'Y -{p-k' - /)2) {p-k' - ly 



{l,p-k'){p-l,k') 



{p-iyj k'^{p-k'Y{{p-k'Y -{p-k' -lY) {p-k'-iY 



{p-k'Y 



1 



7r2 {X + Yf 



167r3 1 3 X2r2 

el 



Y2 



IK 



X2Y2 



du 



Jo 



In '1 



'0 u- 
Inii 



lnii 



X2 



'0 u + 



In u 



{X,A)-iK 



■u - 



A2 



(y,A)+iK 



c.c. 



u - 



A2 



, x2y2 (x,y) + m 



(149) 
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where k = ^JX^Y^ - {X,YY, and 



2{p-k',p-l-k'){p,k') ^^^ {p-k'Y 
p^k'^{p-k'Y{p-k' ^lY k'^ 



(150) 



2(p -k',p-l- k'){p - I, k') {p-k' ~ If 



327r3X2r2 



X- In^ fj + In^ ^ + 2(X, Y) In ^ In 



{p-iYk'\p~k'Y{p~k' - If 



167r3X2r2 



du 



(A,X)+iK 
U- ^ U- 



, . — ^ c.c. In — — , — In ^ 



2G,fc')P" 



In 



2(p - I, k')p- 



k'^[{p-k/f-{p-k' -If] {p-l-k'f k'\p - k' - If - k'"^] 



X In 



{p-l-k'f 



2{l,k'fp-lf 



In 



[p-k'f 



2{p,k'){p-l)-^ ^ Jp-k'f 



k'^[{p-k'f -{p-k' -If] [p-l-k'f k'^[{p - k'f - k' 



1 / 1 X2 I Y'^\ X^ 1 



7r2 1 

6X2 + f3r^ " X2 Jo 



du\a.u 

I 



y2 



u - 



dwlnu 



(151) 



ip-k'-l.k') , (p-k'f 
In 



(p-fc',fc') 



/\2 



p^k'^{p-l- k'f {p-lfk'^{p-k'y 



In 



(P - ^ - fcQ 



1 X2 r 1 , X2 
In 



327r3 r2 



y2 A2 x^ A2 



(152) 



4/p2 



In 



(p-fc') 



n2 



41n(j5-/-A:')Vfc' 



(p-fc')'-(p-fc'-O' (p-l-k'f p2[{p - I - k'f - k'^] 
(p-k'f 41n(p- fc')Vfc'^ 



■In 



1 , ^2 
In 



(153) 



■In 



Y' 



(p-lf[(p-k'f-(p-k'-lf] (p-l-k'f (p-lf[(p-k'f-k'^] 47r3r2 A2 47r3X2 A2 



2(p,p-l-k') ^^ (p-k'f ^ 2(p-l,p-k') ^^^(p-l-k'f 



167r3X2r2 



2 / du 



p^(p-l - k'fk'^ 
In u In 7i 



C.C. 



A2 



(p-k'f k'\p- If 
, x2y2 (x,r) + iK" 

In — — ^ — In - 



(p-lf 



A4 



(X, Y) - in 



iX.Y) ^^2 



(154) 
X2 



167r3X2r2 y2 



Adding the integrals (1481 - (1521 we obtain 



ct^p dH e^(f'^)+^(''^)F2(p,/) 

1 , X2i (X + Yf 



1 

47r3 



1 y2 

^2 A2 y2 A2 



167r3X2y2 j 



if 



du 



In 1 



(A,jf)+»K (A,y)-»K 



327r3X2y2 A2 A2 ^ 87r3X2y2 y2 
Inw 1 1 X2y2 (X,Y)+iK 

-2^"^^" (x,y)-z. 



and therefore 

V^p cf 2; g^(p,A)+^(^,y) [j^reg^^^ ^ ^^^^^ 



(155) 



1 (X,Y) 
87r3 X2y2 

IK 



11 X2y2 



67 7r2 



1 A2 x2 y2 

In — In ■ 



11 



327r3X2y2 A2 A2 487r3 



1 y2 1 X2 
]^l^A^ + 72l^A^ 



167r3X2y2 



du 



In I 



(A,X)+»K (A,y)-»K 

^2 "I ^2 



1 x2y2 (x,y) + iK 
-2^"^^"(x,y)-z. 



(156) 
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Note that the r.h.s. of this equation is finite a,s X ^ Y (taken separately, the contributions of Fi and F2 are singular 
in this limit): 



1 1 
8^X2 



11 .2 67 



(157) 



Using Eqs. (156 1 and (1571 we obtain 



d^p dH e''-P'^^+'<^''^^[Fl'%pJ) + F2{pM (e'^^-'^) - e^(P.^))(e-'(''-''^) - e-'(P-'^^)) 



1 

8^X2^2 

IK 



11 X^Y' 



67 7r2 



1 A2 x2 y2 

In —77 In ■ 



167r3 X2y2 A2 A2 



87r3X2y2 







In 1 



In 1 



(A,X)+iK 



A2 



(A,y)-iK 



A2 



1 X2y2 {X,Y) + iK 



(158) 



Now we turn our attention to last two terms in Eq. (136). Using Fourier transformation 



h 



2i 



In 



(fci + fc2)2fc2 kl 



87r2 



(a;i,a;i2) 

2 ^12i 











/ IK12 





In u 



(3:1,3:12) — 



(159) 



1 (3:2,3^12) + ^^12 

2 ^^2 (a;2,a;i2) - iKi2 



«a;i2i 



167r2a;2 



12 



In^ln^ 

'^12 ^2 



(where 2:12 = xi — X2 and K12 = \/ x\x2 — (xi,a;2)2) one easily obtains 
d^pd^ld^k'ie 



gi(p-fe',x)+i(fe',y) _ gi(p-fe',y)+i(fc',x)-j 



(fc',p-fc') 



{p^iYk''^{p-i-k'Y 



In 



(P - ^ - fc') 



M2 



167r3X2y2 

IK, 

167r3X2r2 



du 



— du 



In u 



Inw 



{X.Y)-iK ' 

In u 

A2 



{X,Y)~iK 



c.c. 



y2 

In u 

(A,y)-iK 

A2 



1 X2 [(A,X) + i/.][(A,y) + »«] 

2 ""^2 [(A,X)~iK][{A,Y)^in] 



1 x2r2 (x,y) + m 

In — — In 7 ^ 

2 A4 (X,y)-iK 



327r3X2y2 y2 



(160) 



Similarly, 
j d'^kid'^k2 e 

167r2 



(fcl + fc2)2fc?fc2 fc2 
1 



i ( (a;i,a;i2) 

a;ii 2 ^12 

)7r ^ X 1 o 



167r2 



167r2 



a;2i 



\Xl2i 



{X2,X12) 



-Xl2, 



•^12 
(2:1,2:12) 



-Xi, 



1 
1 

IK 
1 

IK 



du 



du 



du 



In u 



U T 



In u 



{x2,Xi2)-iK 

u H zn 



Ini 



(xi ,Xi2) — i« 



1 xj {X2,xi2) + in 

- in ^5- In r 

2 x{2 {X2tXi2) - IK 



1 xl {XI,X12) +iK 

- In ^5- In r 

2 a;^2 (2:1, 2:12) - 

1 x\ {xi,X2) + in 

- In In r 

2 XI2 (2:1, X2) — in 



(161) 



327r2a;?2 



.2^2 2 
In ^ In ^ 

-12 ^2 



327r2x2 



2 2 

In^ln^ 
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and therefore 



£f2p£f2;^2^/(-g-»(p-i,X) ^-i{p-l,Y)-^^^i(p-k',X)+i{k',Y) _ gi(p-/c',r)+j(fc\X)^ 

2(p -k',p-l- k'){k',p -l-k') {p-l- k'f 



{p-i)^{p-k')^k'^{p-i - k'Y 



In- 



dui -2 



In u 



In u 



{A,X) + iK " , (A,Y)-tK 



c.c. + In — In ^ ^ 



(X,Y) , X^A^, X2 (X,Y) , y2 1/1 1 \X^Y^ 

— - — - — In In — - — - — In In I 1 In In — 

327r3X2y2 y4 ^2 32tt3X'^Y^ 32tt3\X'^ YV Y^ 



Adding the equations (160) and (|162|) we obtain 



{k',p - k'){p - k')^ - 2{p - k',p- I - k'){k',p -l-k') ^ {p-l- k'f 



(162) 



(163) 



IK 



IGtt^X^Y^ 



du 



Inw 



(p-/)2(p-/c')2fc'2(p-Z-fc')^ 

In u 



In- 



{A,X)+tK (A,y)-»K 



1 x2y2 (x,y) + m 

c.c. 1 in — — — in 



2 A4 {X,Y)-iK 
It is easy to see that the contribution of the last term in Eq. ( |136[ ) is equal to ( |163[ ) so we get 

87r2 



A2 x2, y2 

- In —-7 In 



327r3X2y2 ^2 ^2 



(Tr{[/.;7t})Fig[6].-..- = ^y'd2z[^Tr{[/.[/t}Tr{L/,C/t}-iTr{C/.C/t}] 



X 



+ 2 / £f2p^2;^2^,(g-Kp-i,x) _ g-.(p-i.y))(g,(p-fc',x)+,(fc'.y) _ g«(p-fe',y)+.(fe',x)) 

(k',p- k')(p - k'f -2(p- k',p-l - k')(k',p- l-k') ^ (p-l- k')^ 

In ■ 



{p-l)^{p-k')^k'^{p-l - fc')2 



k'^ 



87r3 



d^z [Tr{U,Ul}TT{UM} - ^Tr{[/,f/t}] 



A2 



11, X2y2 , 67 ^2 



3^" A2 



9 3 



(164) 



Note that the dilogarithms and products of logarithms have canceled. The simplicity of the final result indicates that 
there should be a less tedious derivation but we were not able to find it. 



X. APPENDIX B: CUTOFF DEPENDENCE OF THE NLO KERNEL. 



We will repeat the procedure from Sect. ( Ill C[ ), this time using the cutoff by the slope. 



{KnLO^T^{UxU^}) Shockwave — ^ (Tr{C/a;?7j{})shockwavc — (^LoTr{J7^C/jJ})shockwave 



(165) 



Instead of Eq. ( 18 1 we get 

/•oo /"O 

du dv{A'^{un + x±)A'',{vn + y±)) 

Jo J -oo 

1 2^2 



(166) 



2^ 4 



[j-anlj-bn'l' J (f ad aid f^d P' d Pid p[d P'^ J d^zd^z' j d^ q^d^ q2d^ k^d^ k2 e^(9i+92.^)^-«(fei+fe2,a)i 

4ai(a - ai)t/;"'t/5e-''to-fci'")^-»fe-^^'"')^d.A(QPi + PP2 + qi± + ki^)dy,{api + P'p2 + g2± + ^2^) 
(/? - /3i - /32 + ie)(/3' - /?; - + ie){f3 + - ie){p' + ^a' - ie)[al3s - {qi + q2)l + ^e][a/3's - (fci + fcz)! + ie] 

df^^aiPi + P1P2 + qi±) d^tj.'{aiPi + l^'iP^ + fci±) d^„((a - ai)pi + /32P2 + g2±) dV((a - ai)pi + /3^p2 + ^2^) 
ai/3is — 52^ + ie aiP'^s — k\^ + ie (a — ai)P2S — q2± + — q:i)/?2S — fc|j^ + ie 

r'''''^(Q;pi + qij_, (a - Q;i)pi + q2±, -api - qi± - q2±) T^"" ^ {api + ki±, (a - ai)pi + k2±, -api - kij_ - k2±) 
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where ^ = e~^''i. In this formula comes from the integration over u parameter in the l.h.s. and 

from the integration over v parameter. 

Taking residues at /? = —^a and (3' — — ^a' and P2 — P'2 — —Pi we obtain 

00 /-O 

du I dv{Al{un + XA_)A\{vn + y±_)) 

J —00 



(167) 



fanl fbn'l' 

2^ 4 



aija - ai) jjnn' rjW -i{qi~ki)z^i(q2-k2)z' (gl-L + g2J-)A (fel_L + k21_)x' 



ai/3is — + q:i/3Js — fcjj_ + ie —(a — ai)j3is — q^_^ + —(a — ai)(3[s — k'^_^ + ie 

r^''^(aipi + 9i_L, (a - + q2±,-api - qi± - q2±) T'' " ^ {aipi + ki±, (a - ai)pi + k2±, ~api - ki± - k2±) 



which leads to (cf. Eq. (23 1 



d 
drj 



{Tr{U,Ul})F,,^ = ^Tr{eu,t'ul}n'f^'''j d^ zd' z' U^' U'J,' 



(168) 



"da 



dUo a Jo 



du uu I qid^ q2d^ kid^ k2 



^l(qi,X)^+i(q2,X')^-i(kiX)^-i(k2.Y')^ 



[{qi + q-iY + £.a^] [{ki + fca)^ + ^a^]{qlu + q^u){kfu + k^u) 



{ql - qD^ij — qii{qi + q2)j + -{qi + q2)iq2j {kj - ^2)% — fcii(fci + k2)j + -{ki + k2)ih 



'2j 



9 

47r2 



Tr{eU^t^ Ul }r'^ /''"''' j d^ zd^z'Ur' uf, 



oc /'a 



X ^-^ I da I da' / d^qid^q2d^kid^k2 



{a — Q;')Q;'e*(''i'^~^)^+*('?2:2:-z')i-i(fei,y-2;)j_-i(fc2,J/-z')j_ 

[{qi + qiY + e«2][(fci + k2Y + ia-^]{ql{a - a') + qla'){kl{a - a') + kla') 



S ' ' '2i 2 ir(5'' 2 2 

— {ql - ql) — iqiiiqi + 92)^ h Aqi + q2)iq2j —{kl - kj) — jkuiki + k2)j h -{ki + ^2)4^2^ 

a a a — a JLa a a — a 

(recall that = —2^^). The contribution which is sensitive to the subtraction of (LO)^ is 



drj 



{Tr{U,Ul})^,^^ = ^TT{t-Uj'ul}n'f'"'''' J d^zd'z'Ur'u'J' 



X ^— da da / d qid q2d kid k2 



^i{qi,X)^+i{q2,X')^-i{kt,Y)_L-i(k2,Y')s_ 



'dUo Jo 
a a ~\ 

— ((ji,fci)H 7(92, fe) (91 + 92,^1 + ^2) 

a a — a J 



[{qi + 92)2 + ^a^][{ki + k2Y + ^a^]{ql{a - a') + qla'){kl{a - a') + kla') 

(169) 



The "+" -prescription (47 1 leads to the subtraction 
d 



= ^^^^^,Tr{t^Uj''Ul}n'f"'''J d'zd'z'Ur'uil'J^^^ J d\id\2i'kid^k2 (91 + 92, fci + fc2) 

gi(9i,X)i+i(g2,X')_L-i(*:i,y)i-i(fe2,y')i f f"' da' \ 



a (91,^1) 



(9i,fci) 



[(9i + 92)2+^a2][(fci + fc2)2 + e«'] [io a' ^{ql{a - a') + qla'][kl{a - a') + kla'] qfkf 

"^(92,^2) (92,^2)' 



" da' 
a — a' 



(170) 



Lt-u L[g2(Q,_Q,')+g2^'][fc2(Q,_Q,/)_^/,2Q,/] ^p2 

The details of the upper cutoff in a do not matter since they correspond to changes in the impact factor which do 
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not affect the evolution. For example, 

^ d f°° da 1 
-2C- 



'da' r (gi,fci) 



1 



^ du 



d f da 



da' 



fci) 



{<li,ki) 



dajQ a {qi+ q2Y{ki + k2Y Jq a' l{q\{a - a') + qla'][kl{a - a') + k^a'] qfkf 



(171) 



where the last line is exactly our "rigid cutoff" with "+" subtraction ( 47 1 . 

On the contrary, the details of the upper cutoff in a' are essential for the evolution equation ( 165 1. The contribution 
to {Ki^oTr{IJxU^}) corresponding to the "slope" cutoff ml has the form 



(i^LoTr{;7.{7t})^'°J|§, 



g 

27r2 



Tr{rj7^t''t/]}r"'/''"'''y d^zd^z'Ul'"'uil' J q^d"^ q2d^ k^d"^ k2 



X e 



i(qi,X) ^+i(q2,X') ^~i(kuY) ^-l(k2.Y') ^ 



'da' 



a' 



[qiiki) 



(gi + g2,fci + k2) 
{qi+q2Y{ki + k2Y 
(92,^2) 



{ql + ia'^){kl + ia'^) {ql + ^a''){kj + ^a'^) 



and therefore the difference between the subtractions in "rigid cutoff" ( 170 1 and "slope cutoff" ( 171 1 prescriptions 
can be written as 



^Tr{i'^f/^t''C/]}r"'/''"''' / d^zd^z'Ul'^'uil' / d\id^q2d^kid^k2 e*('?i.-^)^+»fe^^')^-»(fei.^)-L-»('=2,'i")j 



°°da' - 



(qi +92)2(^1 +fc2)2 Jo a' 



{qi,ki 



(92, k2 



{ql + ^a''){kl + {ql + ^a'^m + ^a'^) 



2? 



' da 



{qi + q2,ki + ^2) 



(91,^1) , ('72,^2) 



"da^ 



dUo a [{qi + q2)^+£.a^][{ki + k2)^ + ia^]l qfkf qjkj 

Tr{t''U^t''Ul}r''^f'"''^' J d'^zd'^z'U^'ull' Jd^qid^q2d^kid^k2 e^i<]uX}^+^iq2.x')±~^{kuY)^-^{k2,Y')J 

(92,^2) 



27r2 

^ da' I {qi + 92, fci + ^2) 



a' [ (9l + g2)2(A:i + A:2)2 L(g2+^„/2)(^2^^^,2) (q2+^^,2)(^2+^^,2) 



{qi + 92, ki + /C2) 



(91,^1) , (92,^2) 



[{qi+q2V + ^a'^][{ki + k2r + Ca'^]l qfkj qjk^ 
It is instructive to rewrite this result in Schwinger's notations 
(i^LoTr{C/,?7t}"g'd _ /^LoTr{;7,?7t}^'°P^)Fig|3|. 



27r2 



'da' 



a' 



p2 + ^a' 



{x\g,\z)ur iz\g,\y){z\- 



p2 + ^a' 



jji. 



p^ 



p2 + ^a' p2 ^ ^Q, 



/2 



We see now that the difference between the two regularizations of the longitudinal divergence is given by the difference 
of (LO)^ contributions with cutoffs in a determined by the momenta on the first and on the second step of (LO)^ 
evolution. 
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It is easy to see that for the sum of all diagrams this yields (see eq. ( 55 1 ) 
= all (fzdPz' [Tr{C/,C/t}Tr{C/,C/],}Tr{C/,-[/t} - Tr{C/,[/t[/,,C/tt/,C/],} + [z ^ z')] 



I— , 



\z)-{y\^^\z) (ze|z')-(2/e|z') 



Pi I ,x 



L + 1 + r J L 

Using the integral 



- (z ^ 






z' 






2 


2 , J- 


Pt, 


1 2 





(x\ 

{X,Y) 
47r2X2r2 



1^) 



ln4- 



47r2 ^0 
X2lny2 _y2i^j;^2 

X2 -y2 



(X,F) /-i^^ r(e)r(2-e)(4wzi) 



(X2M + y2u)2-e 

lnX2y2'j +0(e) 



p2 +t' 



(i^LoTr{;7,J7n"g'd - i^LoTr{[/,;7]}^i°P'= 



(fz(fz' 



1 



A2 
Y 



72 {(^'^ 



x2 + y2 X 
X2 -y2 ^72 



(z-z')2X2y2 
2 \-2 

+ 2 



(z-z')2 -y^ ^ 



(172) 



[Tr{;7,C/t}Tr{C/,!7],}Tr{;7,.C/t} - Tr (y^t/.C/],} + z ^ z'] 



(A,y) InX^ - (A,X) Iny^l 

- (y, y') ln(z - z')^ + (y, z - z') in y'^j + a; ^ y 



The NLO kernel for the evolution of color dipoles with respect to the slope is the sum of Eq. (|5| and the correction 
(1721. Note that the correction term (1721 is not conformally invariant (cf. Ref. [24J). This is hardly surprising since 



the non- light-like Wilson line turns into a circle under the inversion x^/x 
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